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COHERENT  SOB-APERTURE  PROCESSING 
TECHNIQUES  FOR  SYNTHETIC  APERTURE  RADAR 


by 

K.H.  Wu  and  M.R.  Vant 


ABSTRACT 


Coherent  sub-aperture  techniques  for  processing  syn¬ 
thetic  aperture  radar  (SAR)  signals  are  described.  The 
techniques  involve  partitioning  the  full  aperture  data,  in 
either  the  time  or  the  frequency  domain.  Each  partition 
or  sub-aperture  is  then  processed  independently  by  the  ap¬ 
plication  of  a  conventional  matched  filter,  or  the  equiva¬ 
lent.  The  low  resolution  images  formed  by  this  process 
are  then  coherently  recombined  to  form  a  full  resolution 
image.  Such  a  processing  scheme  has  the  advantage  of  high 
computational  efficiency  and  provides  the  capability  for 
the  compensation  of  low  frequency  spurious  motions,  even 
when  the  sar  is  operating  in  a  highly  squinted,  high 
resolution  imaging  situation. 

A  detailed  description  of  tvro  coherent  sub- aperture 
techniques  namely,  the  multi-look  matched  filtering 
approach  and  the  step  transform,  are  presented.  The  tech¬ 
niques  are  characterized  and  the  normal  computational 
requirements  are  evaluated.  Computer  simulations,  which 
were  performed  to  verify  the  feasibility  of  the  processing 
schemes,  are  also  described. 


1  .  INTRODUCTION 

In  this  report  two  coherent  3ub-aperture  schemes,  for  compressing 
synthetic  aperture  radar  (SAR)  azimuth  data  are  described.  One  scheme  is 
based  or.  matched  ^tlter^ng  and  <-he  other  on  the  s‘cp  traes'orm.  T r:  pro¬ 
cessing  involves  breaking  up  an  aperture  data  set  into  pieces,  in  either 
the  spatial  or  the  frequency  domain.  Each  piece,  or  sub-aperture,  which 
has  a  fraction  of  the  original  bandwidth,  is  processed  individually  and 
then  the  processed  piece?  are  recombined  coherently  to  restore  the  orig¬ 
inal  bandwidth  and  hence  the  original  resolution.  As  opposed  to  the  usual 
non-coherent  sub- aperture  processing,  where  the  major  purpose  is  to  smooth 
out  speckle  noise  at  the  expense  of  resolution,  the  technique  described 
herein  is  aimed  at  improving  the  processing  efficiency,  within  the  con¬ 
straint  of  having  to  produce  high  resolution  imagery  at  high  squint 
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angles.  The  additional  requirements  that  low  frequency  spurious  target 
motion  must  be  corrected,  and  that  the  whole  system  performs  with  limited 
processor  memory,  are  also  imposed.  It  is  also  assumed  in  the  Zollowing 
discussions  that  the  angle  of  the  antenna,  with  respect  to  the  aircraft  or 
satellite  velocity  vector,  remains  fixed,  i.e.,  the  target  is  not  spot¬ 
lighted.  Most  of  the  techniques  can  also  be  applied  to  spotlight  SARs, 
but  that  problem  is  not  considered  here. 

In  the  first  technique,  matched  filtering  via  fast  convolution,  the 
bandwidth  of  the  recorded  data  is  partitioned  in  the  frequency  domain. 
Each  partition,  or  "sub-aperture",  or  "look",  is  extracted  by  an  FIR 
filter.  (The  terms  "look"  and  "sub- aperture"  are  used  interchangeably). 

Each  look  is  then  compressed  by  a  linear  FM  matched  filter  speci¬ 
ally  tailored  to  overcome  the  complexities  associated  with  large  squint 
angle  and  high  resolution.  In  synthesizing  each  sub-aperture  matched 
filter,  a  piecewise  approximation  of  its  exact  phase  characteristic  at  the 
centre  of  the  look  is  used.  It  is  assumed  that  there  is  little  phase 
deviation  over  the  entire  cross-range  extent  of  the  look. 

As  described  later  in  the  report,  the  processing  requires  that  each 
sub- aperture  data  be  shifted  in  range  to  correct  for  the  effects  of  range 
curvaturs .  The  shifting,  which  must  be  done  by  means  of  an  interpolation 
operation,  requires  that  a  certain  minimum  number  of  rows  or  range  cells 
be  in  memory  simultaneously.  Since  each  3ub-apertuxe  is  shorter  than  the 
full  aperture  it  was  extracted  from,  it  is  frequently  possible  to  meet  the 
minimum  number  of  rows  requirement  with  sub-aperture  processing  whereas 
with  full  aperture  processing  it  is  not.  Even  in  situations  where  the 
minimum  number  of  rows  requirement  is  met  with  the  full  aperture,  it  may 
be  advantageous  to  use  sub-apertures  to  increase  the  interpolation 
efficiency. 

In  situations  in  which  target  motion  must  be  compensated,  the  sub- 
aperture  approach  allows,  at  least  in  theory,  for  piecewise  compensation 
of  the  target  motion  over  the  number  of  sub-apertures.  This  of  course 
assumes  a  suitable  motion  spectrum  for  the  target. 

The  second  technique,  known  as  the  step  transform  [6-9],  is  a 
modified  form  of  the  deramping  technique,  wherein  small  deramp  references, 
whose  lengths  are  a  fraction  of  the  full  aperture,  are  used  to  extract 
sub-apertures.  Upon  deramping  each  look,  the  data  set  becomes  an  ensemble 
of  CW  signals  whose  frequencies  correspond  to  the  relative  spatial 
locations  of  the  targets.  Then,  the  target  spectra  can  be  extracted  using 
a  filter  bank,  which  can  be  realized  by  the  discrete  Fourier  transform. 
The  process  is  repeated  for  the  other  sub-apertures,  and  the  extracted 
spectra,  from  the  sub-apertures  that  are  associated  with  a  given  target, 
are  added  coherently  to  form  an  image  of  the  target  at  full  resolution. 
The  major  advantages  of  thi3  technique  include  a  small  requirement  for 
memory  and  a  smaller  signal  loss  than  with  the  conventional  deramping 
technique.  This  report  extends  the  original  concept  of  the  step  transform 
[6-9)  and  shows  how  such  effects  as  range  curvature  and  cubic  phase  error 
can  be  overcome. 
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In  the  following  sections  details  of  the  different  principles 
underlying  the  technique  of  sub-aperture  processing  via  matched  filtering 
is  given.  Then,  the  general  form  of  the  radar  return,  including  such 
effects  ac  range  curvature  and  cubic  phase  error  is  formulated,  followed 
by  a  description  of  the  image  quality  degradations  that  result  if  these 
effects  are  left  uncompensated.  Next,  the  general  form  of  the  sub¬ 
aperture  (multi-look)  matched  filter  is  synthesized,  using  a  piecewise 
approximation.  This  is  followed  by  an  analysis  of  the  projected  per¬ 
formance  of  the  technique,  based  on  the  computer  simulation.  Next,  an 
overview  of  the  various  deramping  techniques  is  given  and  finally,  the 
computational  requirements  for  the  approaches  are  calculated,  and  their 
efficiencies  compared. 

2 .  BASIC  CONCEPT  OF  COHERENT  SOB-APERTORE  PROCESSING  VIA  MATCHED 

FH.TERIHG 

The  principle  of  sub-aperture  processing  is  based  on  the  fact  that, 
provided  the  frequency  variation  with  respect  to  time  is  monotonic,  and 
the  amplitudes  in  the  data  change  gradually,  there  is  a  one-to-one  mapping 
between  the  time  and  the  frequency  domains.  This  phenomenon  is  known  as 
the  principle  of  stationary  phase  [ 3 ] .  How  this  principle  relates  to  SAR 
processing  is  explained  next.  Point  targets  paraded  on  a  straight  line 
parallel  to  the  radar  flight  path  will  give  exactly  tne  same  phase 
history,  and  will  differ  only  in  their  respective  time  delays.  If  the 
collected  target  return  ensemble  is  Fourier  transformed,  the  frequency 
responses  of  all  targets  will  be  superimposed  on  the  same  frequency  range, 
but  within  the  return  from  each  target  there  will  be  embedded  a  linear 
phase  tern.  ( in  addition  to  other  inherent  higher  order  phase  terms ) ,  which 
signifies  its  spatial  location.  A  look  extraction  process  in  the 
frequency  domain  is  equivalent  to  simultaneous  lock  extractions  for  all 
targets  illuminated  by  the  same  portion  of  the  antenna  pattern.  The 
portion  of  the  spatial  look  thus  extracted  depends  on  the  frequency-space 
relationship,  which  in  turn  is  governed  by  the  geometrical  complexity  of 
the  imaging  scenario.  The  relationship  between  frequency  and  space  is 
greatly  simplified,  indeed  it  is  one-to-one,  if  the  principle  of 
stationary  phase  holds. 

In  order  to  facilitate  the  illustration  of  the  concept,  we  assume 
that  the  Doppler  phase  history  traced  by  a  point  target  is  purely  quadrat¬ 
ic  (i.e.  linear  FM)  over  a  full  aperture,  of  duration  T.  This  typifies 
the  side-looking,  low  resolution  case.  The  form  of  this  one-dimensional 
signal  is 

jTTBt2 

u (t )  »  rect ]  e  ,  *  1 ^ 

T 

i-  t  1 

where  B  is  the  bandwidth  in  Hz,  and  rect  [—  I  is  a  rectangular  function  of 

T  ' 

duration  T.  Assuming  the  time-bandwidth  product  is  sufficiently  large, 
the  Fourier  transform  of  u(t)  can  be  approximated,  using  the  principle  of 
stationary  phase  [3]  as: 


The  rect[  ]  can  be  thought  of  as  an  ideal  low  pass  extraction 

2q+l 

filter  with  bandwidth  B/NL  and  centred  at  — j— B . 

L 

If  the  qth  loox  spectrum  is  shifted  to  baseband  by  substitution 
f  +  f '  +  2H11b,  (4) 

2NL 

it  becomes , 

Uq(f)  “  DqU  +  B)' 

L 

NLf  B  2NL 

*  rect[ - j  e  , 

B 

-£H(lS±ls)r  .J2H(i3lL)2B2 

rNLf  •  .  B  B  2NL  B  2Nl 

-  rect[ - J  e  e  e 

B 
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they  must  be  taken  out  before  the  universal  matched  filter  is  applied, 
when  this  is  done, 

+  j££(1SIL)2b2 

U"q(f)  -  0'q(f')  e  B  2nL  b  2nL  , 

-jllf  -2 

rNtf*  B 

■  rectf - ]  e  ,  (6) 

B 

is  obtained.  The  matched  filter  for  C(f')  is 

<3 

J^If'2 

M(f')  »  rect[Nkf  ]  e  0  < 7 ) 

B 

After  multiplication  by  the  matched  filter  in  the  frequency  domain,  the 
compressed  signal  is 


Rq(f')  -  U-(f')  M(f')  -  rect[^l£-]. 


After  frequency  shifted  back  to  its  orignal  frequency  position,  using  (4) 
for  the  relationship  between  f  and  f,  the  correlated  signal  becomes 


(9) 
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The  same  process  can  be  applied  to  all  the  other  looks,  and  then  the  looks 
cam  be  summed  or  concatenated  together  to  form  the  original  bandwidth 
signal, 

N- 

—  -1 
2 


R  ( 


M.  B  2Nr. 


q»- 


.  r  f ' 

rect |— 
B 


no ) 


The  compressed  pulse,  or  the  point  spread  function  of  the  point  target,  is 
obtained  by  taking  the  inverse  Fourier  transform  of  R(f'),  i.e„. 


r(t)  ■  B  sinc(Bt ) 


(11) 


Alternatively,  the  summation  can  be  performed  in  the  time  domain 

r- 

r(t)  -  l  F_1{rect[— 

Nj  B  --L 


-1  12*  (  )Bt 

F  sincffi  - 


NL  nl 


(?) 


B  sinc(Bt)  • 


(12) 


The  above  analysis  is  illustrated  in  Figure  1 ,  where  two  point 
targets  located  at  t®0  and  t^t1  and  at  the  same  3lant  range  are  assumed  to 
be  present.  A  Fourier  transform  over  the  entire  time  domain  results  in 
superimposed  spectra  of.  the  two  targets,  since  they  both  have  the  same 
phase  characteristic.  Moreover,  there  is  a  one-to-one  correspondence 
between  the  time  and  the  frequency  domains,  as  illustrated  by  the  mapping 
of  two  corresponding  strips  of  the  time  signals  into  a  single  superimposed 
strip  in  the  frequency  domain.  The  spectrum  is  then  partitioned  into 
equal  looks.  Each  look  is  individually  processed.  All  the  processed  look 
spectra  are  then  concatenated  and  inverse  Fourier  transformed  to  form  the 
final  compressed  pulses.  Since  the  system  is  linear,  an  ensemble  of 
target  returns  can  be  processed  simultaneously.  The  processing  procedure 
is  summarized  in  Figure  2,  where  the  coherent  look  summation  takes  place 
in  the  time  domain. 


signal  amplitude. 


J - - — i - -  time 

le 


Pig.  2  -  Coherent  sub-apertare  processing  procedure  for  an  ideal 
linear  PM  signal. 


Other  than  the  fact  that  the  foregoing  approach  can  ease  the 
computer  memory  requirement  in  some  cases,  it  is  obvious  that  sucn  a 
processing  scheme  would  be  redundant  if  the  input  signal  is  as  simple  as 
that  given  by  (1).  However,  complications  arise  under  different  operating 
conditions  such  as:  high  squint  angles,  long  integration  periods,  and 

spurious  target  motions.  In  order  to  combat  these  complications,  with 
little  loss  of  processing  efficiency  and  small  increase  in  memory 
requirements,  it  may  be  necessary  to  fine-tune  a  matched  filter  for  each 
segment  (look  or  sub-aperture)  of  the  input  data.  The  inter-relationship 
between  the  sub-aperture  and  this  fine-tuning  procedure  will  be  discussed 
further  in  the  following  sections. 
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3  .  GENERAL  FORM  OP  RECEIVED  SIGNAL 


The  received  signal  in  the  azimuth  dimension,  encompasses  more  com¬ 
plications  than  the  one-dimensional  linear  FM  signal  described  in  the  pre¬ 
vious  section.  when  the  antenna  is  scruinted  away  from  the  side-looking 
position,  abnormalities  in  phase  and  time  of  arrival  of  the  radar  return 
sional  become  a  problem.  In  general,  higher  order  phase  terms  (e.g.  the 
cubic  term)  may  severely  degrade  the  compressed  pulse  width  and  the  inte¬ 
grated  sidelobe  ratio.  If  the  time  of  arrival  of  the  radar  return  signal 
and  the  digitization  timing  are  not  synchronized  then  the  range  compressed 
target  history  associated  with  a  point  target  is  no  longer  confined  to  a 
single  row  of  digitized  azimuth  data,  but  instead  spans  several  range 
cells.  This  phenomenon  is  commonly  known  as  range  migration.  A  thorough, 
treatment  of  these  effects  can  be  found  in  [2]. 


Next,  we  will  describe,  in  terms  of  the  slant  range  and  squint 
angle,  a  general  form  of  the  received  signal.  This  formulation  will  pro¬ 
vide  for  a  better  understanding  of  the  various  complications  and  hint  a- 
ways  to  suppress  image  quality  degradations,  while  still  maintaining  a 
reasonable  degree  of  computational  efficiency. 


With  reference  to  the  flight  geometry  depicted  in  Figure  3,  the 
distance,  r(s),  between  the  radar  and  a  point  target  on  the  ground,  for  a 
flat  Earth  model,  is 


r(s) 


23^003^, 


M3) 


where  s  is  the  distance  along  the  radar  flight  path,  rp  (»r(o))  is  the 


Fig.  3  -  Air bo rue  SAR  geometry  with  fixed  antenna  squint 
angle  (assuming  flat  Earth). 
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slant  range  at  the  centre  of  the  aperture,  and  hp  i3  the  squint  angle, 
measured  between  the  slant  range  and  the  aircraft  velocity  vector,  as 
shown  in  Figure  3, 

For  ease  of  analys  »,  the  right  hand  side  of  (13)  is  expanded  into 
a  Taylor's  series,  i.e. 


where  ag  rp, 


»1 

*2 


cos  Hp, 


sin  He 


2r, 


.  2 

coshpsin  Pp 

2r£ 


and. 


a4 


2  2 
sin  r>p(5coB  nF-l ) 

3 

rF 


ar.d  the  futiu  of  the  received  signal  car.  be  expressed  as  a  two-dimensional 
function,  i.e.. 


.4*  /  » 

-j— r(s) 

gD(t',s)  ■  4>[t'  -  —  r(s)  ]  e  * 

C 

where  C  is  the  velocity  of  propagation. 


(15) 


t'  -  t-mA, 

is  the  time  measured  from  the  start  of  the  mth  pulse,  and  represents  the 
along  range  dimension, 

s  *  mV^A, 

i3  the  distance  travelled  by  the  sub-aircraft  point  along  the  aub-aircraft 
track  during  m  interpulse  periods  of  length  A  sec.,  and  represents  the 
azimuth  (cross-range)  dimension,  VTO  is  the  radar  platform  velocity,  ). 
is  the  transmitted  carrier  wavelength  in  metres,  and  \p ( t )  is  the 
compressed  range  profile.  The  time  delay  embedoed  in  signifies  the 
range  curvature  variation  over  the  aperture.  The  scattering  magnitude  and 
the  antenna  pattern  are  assumed  to  be  unity. 

The  series  given  in  (14)  converqe.  In  most  practical  cases,  when 
calculating  the  amount  of  ranqe  migration  in  y,  only  the  terms  up  to  the 
quadratic  need  to  be  included.  However,  when  calculating  the  phase  term, 
which  must  be  accurate  to  within  0.75u  at  aperture  edge  with  heavily 
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windowed  data,  the  cubic  term  sometimes  has  to  be  included  (see  Appendix 
A).  With  these  modifications  the  general  form  of  the  received  signal 
becomes 


2  i|)[t'  -  —  (a0  *  a^s  +  a 2s  ^ ) 


,4n , 

_  j— ”[  ao  +  als 

e  * 


7  i 

a2S  +  ajs  ] 


(16) 


The  envelope  term  in  (16)  shows  that  the  locus  of  the  ridge  of  the 

range  compressed  data  is  defined  by  a  parabola  ,t '  ^ — (a  +aiS+a29‘1),  in  the 

C 

2 

(t*,j)  domain.  The  tilted  linear  locus,  due  to  the  a23  term,  is  known  a3 
range  curvature  [9].  In  sguinted  mode  SAR,  the  amount  of  range  walk  can 
span  several  range  cells,  thereby  prohibiting  the  use  of  one-dimensional 
processing.  In  practice,  range  walk  can  be  avoided  by  acquiring  data 
along  the  line  a^a^s,  this  simplifies  the  argument  of  i|/;  and  by 
demodulating  the  data  in  azimuth  to  remove  the  a-js  phase  term.  Thus 
gT(t',s)  can  be  rewritten  as 


grft'.s)  *  *[f  -  i  a2a2] 


— — f  Sq  +  a2s  t  a s  ] 

S  "  .  (17) 


The  acquisition  scheme  is  implemented  by  triggering  the  A/D  conversion 

—  (ju+ais)  sec.  after  an  FM  pulse  is  transmitted.  As  the  radar  plat- 
C 

form  advances  in  the  s-dimension,  the  A/D  triggering  delay  time  also 
changes  accordingly.  In  doing  so,  individual  target  data  are  placed  on 
arectangular  grid  of  a  two-dimensional  computer  memory  array,  and  thereby 
minimizes  range  migration  in  ths  range  dimension.  This  is  shown  in  Figure 
4.  Tha  recorded  data  look  as  if  they  were  taker,  in  the  side-looking 
position  (nj»  «=  90*).  However,  there  is  a  major  complication:  targets 
with  different  true  ranges,  rp,  and  thus  requiring  different  phase 
compensation,  are  mapped  into  the  same  processing  cell.  The  processing 
3cheme  outlined  here  does  not  allow  different  compensation  to  be  applied 
to  the  various  targets  if  they  lie  in  the  same  cell.  Therefore,  the 
region,  over  which  the  same  compensation  can  be  used,  must  be  large  enough 
to  accommodate  the  change  in  true  range  with  azimuth  as  we  move  from  one 
end  of  the  cell  to  the  other. 

3  . 1  Spatlal-to-Spectral  Mapping  in  Azl  aith 

Since  the  look  extraction  and  matched  filtering  are  performed  in 
the  frequency  domain,  it  is  necessary  to  know  the  form  of  the  range 
curvature  in  that  domain.  The  mapping  of  range  curvature  between  space 
and  frequency  will  not  alter  the  overall  parabolic  shape  of  the  curvature 
significantly.  Based  on  the  principle  of  stationary  phase,  it  can  be  seen 


11 


DATA  IN  MEMORY  IF 
ACQUIRED  parallel 
TO  a,s 


Fig.  4  -  Two  data  aoquistion  scheau,  (1)  parallel  to  flight  path, 
(it)  parallel  to  ajs.  (a)  Show*  flight  grot  tries  for 
three  targes ,  A,  B,  C.  (b)  Shows  the  corresponding  data 
■naory  arrangement  for  the  two  schemes.. 


that  the  most  important  term  in  the  mapping  is  the  quadratic  phase  term, 
4ir  2 

-r —  a2s  .  This  reinforces  our  notion  that  the  shape  is  parabolic.  However, 
^  4  IT  ^ 

it  is  not  that  simple.  The  cubic  phase  term,  y-  a3s',  introduces  a  non¬ 
linearity  into  the  mapping,  which  slightly  distorts  the  parabolic  shape  of 
the  range  curvature.  The  form  of  the  mapping  i3  discussed  in  the  next 
paragraph. 


The  point-to-point  correspondence  between  the  two  domains  is 
obtained  by  differentiating  the  phase  term,  8(s),  given  in  (17): 


0(3) 


47T 

X 


l*o  +  ®2S  +  a3s 


(18) 


so 


i  de(s) 


2_ 

X 


(2a2s 


+ 


3a  33 


ZT 

;  * 


f 


f  (s) 


2ir  ds 


(19) 


E 


0 


The  solution  to  this  quadratic  equation  in  s  is 

i  -2a?  +  /  4a?2  -  6a3Xf  , 

|  - L - a - i —  ,  for  a  3  *  0; 

633 


(21  ) 


Xf 

4a2 


for  a  3  *  0. 


The  f-s  curve  given  by  (19)  is  sketched  in  Figure  5, 


Fig.  5  - 


d  1  ( f )  ■ 


Dependence  of  frequency  f  on  distance  s.  Look 
partitions,  denoted  by  LK-2,  Ut-'l,  LKO,  UC1 . 

r 

282  r-2a2  ♦  /  **22_zJ*\x± ]  ,  tor  a 3  *  0  ? 

CT  6a3 

I  (22) 

\Z  2 

- fz  ,  for  aj  ■  0. 

8a2CT 


where  T  is  the  range  sampling  interval.  In  general,  a3  is  small, 
can  be  approximated  as 


and  d 
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l2  2 

<3 1  (  f )  =«  - f  . 

S&ZCT 

which  is  a  parabola,  as  noted  before. 


(23) 


4.  ONE-DIMENSIONAL  AZIMUTH  PULSE  COMPRESSION  WITHOUT  RANGE 

CURVATURE  COMPENSATION 

The  exact  matched  filter  for  the  SAP  given  by  (17)  is  two- 
dimensional  [2 |.  Two-dimensional  data  processing  is  characterized  by  its 
massive  data  size,  and  its  efficiency  is  contrained  by  the  available 
computer  memory  space.  Therefore  it  would  be  advantageous  if  one¬ 
dimensional  processing  could  be  used  instead.  This  can  be  achieved  under 
either  of  two  circumstances: 

(i)  if  the  range  and  azimuth  signals  can  be  decoupled  by  break¬ 
ing  up  the  aperture  into  small  pieces,  each  one  exhibiting 
negligible  curvature, 

or 

(ii)  if  the  range  curvature  is  small  enough  to  be  disregarded. 

In  this  section,  we  will  assess  the  degradation  due  to  uncompensated  range 
curvature  and  in  Section  5  we  will  consider  the  decoupling  approach. 

The  dominant  deterioration  caused  by  range  curvature  is  the  broad¬ 
ening  of  the  mainlobe  of  both  the  azimuth  and  the  range  compressed  pulse. 
We  will  attempt  to  characterize  the  mainlobe  broadening  effects  in  range 
and  azimuth  in  a  simple  manner,  while  still  keeping  the  description 
general.  Because  of  the  non-linear  and  three-dimensional  nature  of  the 
scenario,  namely,  range,  azimuth  and  signal  amplitude,  a  more  precise 
characterization  can  only  be  obtained  by  empirical  means,  and  then  only 
for  a  limited  set  of  parameters. 

4 . 1  Mainlobe  Broadening  in  the  Azimutn  Dimension 

A  typical  range  compressed  profile  exhibiting  range  curvature  is 
shown  in  Figure  6.  Profiles  of  the  uncompressed  azimuth  signal  in  the 
spatial  and  the  frequency  domains  are  shown  in  Figures  7(a)-(c).  The 
envelopes  of  these  profiles  are  shaped  by  coupling  between  the  range  and 
azimuth  signals.  Antenna  pattern  effects  have  been  neglected.  Provided 
that  a  Hamming  (or  similar)  window  is  used  during  range  compression  and 
that  the  range  curvature  is  approximately  parabolic,  the  overall  3hape  of 
the  profiles  will  be  insensitive  to  practical  parameter  variations. 

In  this  characterization,  the  azimutn  profile  passing  through  the 
origin  (i.e.  at  s-0,  and  f-0)  was  found  to  be  the  major  energy  contributor 
(see  Figure  7(a)).  This  centre  profile  was  used  to  examine  the  effects  of 
range  curvature,  A  generalized  parameter  P  was  used  to  characterize  the 
extent  of  the  cuvature.  R  is  the  ratio  of  the  full  bandwidth,  Np  ceils, 
measured  along  the  range  curvature  locus,  to  the  -vdB  (half  amplitude) 


14 


<0! 


Pig.  6  -  Signal  history  of  a  point  target.  (a)  Range  curvature 
with  oaepreased  range  profile.  (b)  Perspective  view  of 
range  curvature. 


bandwidth,  2n '  cells,  measured  along  the  middle  of  the  range  cell,  i.e.. 


(24) 


The  parameters  R  and  n'  are  shown  in  Figure  8.  Note  that  n' 
locates  the  -6dB  points  in  range  and  azimuth.  The  offset  of  the  range 
profile  is  given  by  d,  which  is  measured  from  the  centre  of  curvature  (see 
Figure  8).  The  parameter  d  can  be  used  to  relate  n'  to  the  other 
parameters.  First,  d  is  expressed  in  terms  of  the  number  of  range  cells 
of  curvature: 


(25) 


where  the  1.81  factor  is  the  -6dB  mainlobe  width  of  a  Hamming  weighted 
pulse,  expressed  in  terms  of  range  resolution  bins  (see  14],  p.55),  Yj; 
is 


Yr 


sampling  frequency 
range  signal  bandwidth 


where  yR>1 


(26) 


RANGE  CELL  No.  1 


n—  *—  I  — . I  ■  - '  T - — T-'— '  ■■  . — —  ■ — *'■■■—  I 

-2000  -1500  -1000  -500  C  500  1000  1500  2000 

SAMPLE  NUMBER  m 
AZIMUTH  SIGNAL 


uj 

a 

Z> 

£ 

Z 

a 

•< 

s 


•2000  -1500  -1000  -500  0  500  1000  1500  2000 


SPECTRAL  NUMBER 
AZIMUTH  SPECTRUM 


(a) 


Fig.  7  -  Slices  of  uncompressed  azimuth  profiles  (compressed  in  the 
range  dimension )  in  space  and  frequency  domains „  In  the 
signal  profiles,  the  horizontal  axes  denote  the  sasples 
collected  at  the  pulse  repetition  frequency.  The  vertical 
axes  denote  the  real  part  of  the  signal.  The  azimuth 
spectra  are  obtained  by  taking  FFT  (4096  points  in  this... 


-2000  -1500  -1000  -500  0  500  1000  1500  2000 


SAMPLE  NUMBER  m 
A2IMUTH  SIGNAL 


-2000  -1500  -1000  -500  0  500  1000  1500  2000 


SPECTRAL  NUMBER 
AZIMUTH  spectrum 
(6) 

Pig.  7  (cont. )  - 

case)  on  the  received  ccaaplez  signal.  The  vertical  axes 
denote  the  aegnitndea  of  the  spectra.  (a)  corresponds  to 
the  profiles  through  the  centre  of  the  range  curvature 
signal  (see  figure  6),  (b)  and  (c)  are  the  profiles  in  the 
next  two  adjacent  range  cells. 
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RANGE  CELL  No.  a 


•2000  -1500  -1000  -500  0  500  1000  1500  2000 


SPECTRAL  NUMBER 

azimuth  spectrum 
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Flq.  17  {cent. )  - 
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Fig.  8  —  Range  (ooapreosed)  -  asiaotb  (nncoqpreseed)  coup  lag . 


d  -  a2n,2ASA2 

-  n '  2  ASA2 , 

2rF 


(27) 


where  ASA  is  the  azimuth  sampling  interval  in  metres.  The  desired 
expression  for  n‘  am  now  be  obtained  by  equating  (25)  and  (27),  i.e.. 


0.905yrASr 


*inS  n’2AS2  , 
2rP 


(28) 


and  solving  for  n': 

n-  -  !i 781  ^rFASB 

/  AS2sin2np 


(29) 


This  expression  for  n'  can  then  be  substituted  in  (24)  to  obtain  the 
desired  expression: 

n_ASa  sinnF 

R  “  - - -  (30) 

2  /TTSl 


The  azimuth  resolution,  as  derived  in  Appendix  B,  is  given  by 
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p  3  - Llifl -  .  (31) 

2NpASAsinnj? 

Thia  can  be  substituted  in  (30)  to  further  simplify  the  expression  for  R, 
i.e. , 

r  =  /_££_  .  (32) 

P  /  YrASr 


The  larger  the  curvature,  the  larger  the  value  R  will  be.  For 
small  curvature  (R<1),  the  -6dB  point  is  beyond  the  end3  of  the  azimuth 

Nd  -  |  Np 

aperture  bounded  by  — ,  i.e.,  n'  <  — — .  However,  even  in  this  situation, 

2  1  1  —  2 

R  remains  a  valid  parameter  for  the  characterization. 


Computer  simulations  of  the  amount  of  mainlobe  broadening  were  per¬ 
formed  for  different  values  of  R.  The  results  are  plotted  in  Figure  9. 
The  percent  broadening  of  the  compressed  mainlobe  width,  measured  at  the 
-1 OdB  points,  was  defined  as 


Fig.  9  -  Scpirical  azimuth  beam  broadening. 
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W-  -  -■U.~wg  X  100%  (33) 

wc 

where  wu  is  the  mainlobe  width  without  range  cur  mature  compensation,  and 
wc  is  the  mainlobe  width  with  range  curvature  compensation.  Before  wu 
and  wc  were  measured,  the  compressed  mainlobes  were  interpolated  by  a 
factor  of  100.  An  empirical  equation  relating  WA  and  R,  was  developed 
to  fit  the  measurements: 

WA  =  3.4R3-1,  for  R  <  2,  WA  <  2.7%.  (34) 

In  the  above  computer  simulations,  a  Hamming  window  was  applied  in 
the  frequency  domain  over  the  uncompensated  range  curvature  azimuth  data 
set  before  the  inverse  Fourier  transform  was  taken.  Therefore,  Np  is  no 
longer  a  measure  of  the  bandwidth  associated  with  a  single  point  target 
but  is  instead  a  measure  of  the  length  of  the  weighted  aperture. 

4.2  Mainlobe  Broadening  in  the  Range  Dimension 

With  significant  range  curvature,  energy  is  spilled  over  onto  the 
concave  side  of  the  range  curvature  thereby  broadening  the  already  com¬ 
pressed  range  profile.  If  a  Hamming  window  is  applied  in  the  frequency 
domain  in  azimuth,  energy  at  the  end  of  the  curvature  i3  highly  atten¬ 
uated,  and  the  spill-over  in  the  range  dimension  is  small.  Computer 
simulations  were  used  to  demonstrate  that  under  nor.nal  operating  con¬ 
ditions  the  rate  of  mainlobe  broadening  in  range  is  slower  than  that  in 
azimuth. 


A  parameter  x,  which  is  the  number  of  range  cell  crossings  at  the 
end  of  the  aperture,  normalized  by  the  -6dB  range  mainlobe  width,  was  used 
as  a  simple  and  effective  way  of  characterizing  the  range  broadening. 
This  parameter  x  is  given  by: 

_  number  of  range  cell  crossings 
x  ■  —  t 

1.8lYR 

a2ASA2  (  )2 

1.81 Yr 


from  (31 ) , 

0.0338A2rp  . 

x  -  - £  .  (35) 

Z 

YrP 

Computer  simulations  of  the  amount  of  range  mainlobe  broadening,  expressed 
as  a  percentage,  were  performed  for  different  values  of  x;  the  results  are 
plotted  in  Figure  10.  An  empirical  equation 


WR  =  1.46x2-562,  for  x  <  2,  wr  <  9% 


(36) 
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was  found  to  fit  the  curve. 


Fig.  10  -  R^irical  range  bean  broadening. 


5.  A  PUSCEWISK  QUADRATIC  APPROXIMATION  -  COHERENT  SUB-APKRTUKK 

PROCESSING 

In  this  section,  a  generalized  coherent  sub-aperture  processing 
technique  is  described.  This  technique  overcomes  the  processing  problems 
caused  by  the  non-linearities  that  arise  whan  imaging  at  a  high  squint 
angle  and/or  to  a  high  resolution.  The  reasons  for  favouring  the  sub¬ 
aperture  approach  over  the  continuous  straightening  approach  [11]  are  two¬ 
fold; 

i)  spurious  low  frequency  motion  effects,  if  they  exist,  can  be 
corrected  over  small  areas  by  applying  different  compensa¬ 
tion,  e„g.  azimuth  FM  rates,  to  each  sub-aperture;  and 

ii)  The  sub-aperture  technique,  with  its  piecewise  correction, 
requires  much  shorter  row-lengths  than  the  continuous 
f-traiqhtening  approach,  thuu  the  acrosa-row  operations  re¬ 
quired  for  the  interpolation  during  range  curvature  correc¬ 
tion  car.  be  nuch  more  easily  accommodated,  i.e.,  mere  rows 
car.  I*  fit  in  memory  simultaneously. 

Ine  major  goal  is  to  devise  a  3ub-aperture  processing  scheme  chat 
ct  '  .'tmevo,  without  resorting  to  two-dimensional  processing,  the  distor¬ 
tions  utributed  to  the  non-linear  effects  previously  described.  As  with 
m  at  problems  involving  non-linearities,  piecewise  approximation  was 
deen »d  to  be  the  sin, pleat  approach.  Each  piece  of  the  signal  was  assumed 
to  be  linear  up  to  a  certain  degree  of  error  tolerance.  Figure  11  illus¬ 
trate  i  a  slant  range  plane  containing  the  flight  path  and  a  point  target 
of  interest.  The  light  path  is  broken  up  into  several  segments  to  signify 
different  sub-apert-res.  Each  sub-aperture  i»  extracted  by  an  FIR  filter 


and  i.3  Chun  demodulated  by  a  matched  filter,  especially  designed  for  it. 
The  extraction  filtering  and  the  matched  filtering  can  be  combined  in  the 
frequency  domain  and  die  two  processes  performed  simultaneously.  All  the 
sub-apertures,  except  the  two  end  ones,  which  may  only  be  partly  filled, 
have  equal  bandwidt'. 


Pig.  11  -  A  slant  plane  containing  flight  path  and  two  point 
targets,  a  and  b.  Note  the  corresponding  look 
designations  (denoted  by  LX )  are  displaced  by  the  sane 
waul  am  the  separation  between  two  targets.  The  look 
geometries  are  congruent. 


In  order  to  define  the  matched  filters  for  each  sub-aperture 
several  other  things  must  first  be  done.  The  look  centres  in  the 
frequency  domain  must  be  located,  and  their  frequency  domain  locations 
translated  into  corresponding  spatial  domain  ones.  This  ties  down  the 
values  of  effective  squint  angle,  hq,  of  the  qth  look,  and  slant  range 
to  the  look  centre,  rq,  which  are  required  to  define  the  matched 
filter.  Next  the  phases,  over  the  look  and  over  the  local  section  of  the 
full  aperture,  are  matched.  The  look  phase  contains  terms  of  up  to 

quadratic  order,  whereas  the  full  aperture  phase  contains  terms  of  up  to 
cubic  order.  The  phase  matching  must  be  done  in  order  to  maintain 

coherency  over  the  phase  jumps  between  the  looks.  Range  curvature  com¬ 
pensation  must  also  be  apnlied;  in  this  case,  as  already  discussed,  the 

curvature  compensation  is  to  be  piecewise.  Finally,  the  demodulation 
whicn  removed  the  ajs  term  during  data  acquisition  must  be  accounted  for. 
After  all  these  operations  are  completed  the  matched  filter  can  be 
defined. 

Referring  back  to  Figure  1 1 ,  each  look  can  be  thought  as  inde¬ 
pendent,  with  its  own  look  centre,  slant  range  and  squint  angle.  The 
locations  of  the  look  centres  are  based  on  the  number  of  looks  required, 
or  equivalently,  the  bandwidth  of  the  extraction  filter.  In  the  present 
scheme,  the  entire  frequency  axis  is  extracted  even  though  part  of  the 
spectrum  has  no  signal. 

In  the  frequency  domain,  the  look  centre  of  a  sub-aperture  is  given 

by. 
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2q-1  NL  N, 

f_  “  -  - *1  <  q  <  -  ,  NL  is  even,  (37) 

q  2NL  2  2 

where  q  denotes  the  look  index,  and  Be  is  the  bandwidth  of  a  pre¬ 
designed  extraction  filter.  3e  is  a  function  of  the  number  of  looks, 
N^,  such  that  BeNL=0.9fg,  where  fs  is  the  sampling  frequency.  A 
bandwidth  of  O.lfs  is  reserved  as  guard-band  near  che  Nyguist  frequency. 


(3!?) 


(40) 


Equipped  with  hq  and  rq,  and  the  previously  derived  quadratic  approxi¬ 
mation,  we  can  synthesize  the  phase  of  the  acquired  data.  Within  the 
q'-!l  sub-aperture ,  the  phase  is, 

Oq(s')  -  -  ^  (a^  +  a,qs'  +  a2qs'2),  (4i) 

where  s'  is  the  spatial  co-ordinate  with  its  origin  located  at  the 
sub-aperture  centre, 

«oq  "  tq,  (42) 


a  i  q  “  -cos  riq 


(43) 
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*2  q 


sin2nc 

2r„ 


Prom  (16), 
9(3) 


the  phase  of  the  original  data  is  given  by, 

4ir  ->  3 

-  -  (»c  +  *13  +  a2s~  +  a3s  ). 


(44) 


(45) 


The  phases  0q(s')  and  0(s)  must  match  very  closely  over  the  look,  i.e., 

8q(s')  »  9(s),  for  ft  *  3q,  s'  "0, 

and  (46) 

0a(s')  3  6(s)  ,  for  3_-  ?s^sa-iis_,  i.e.,  -  ^3.  ZB '  ?  ^_2_ , 

q  q  2dL  q  2NL  2Nl  2Nl 

where  Lg  is  the  full  synthetic-aperture  length.  Since  s'  ■■  s-sq, 
0q( s ' )  can  now  be  written  as. 


0q(B)  -  -  —  [aoq  *  alq(s-sq)  +  a2q(s-sq)2], 


(47) 


Because  amplitude  matching  is  unimportant,  we  can  set  the  amplitude  to  be 
unity  without  piecewise  approximation,  and  rewrite  the  piecewise 
synthesized  phase  signal  ovtr  the  sub-aperture  as 


t *oq  +  aiq^-3q)  +  a2q(s-sq)2] 

e 

(4e ) 


i|>(t)  represents  the  profile  of  the  range  compressed  pulse  over  the  range 

,  2  ?  ’ 

dimension,  and  the  variable  It'-  —  a2s  |  represents  the  continuous  range 

C 

curvature,  as  a  function  of  position  in  the  azimuth  aperture.  As 
described  in  Section  3,  range  curvature  correction  is  sometimes  required 
to  prevent  range  and  azimuth  broadening  from  occurring. 


It  range  curvature  compensation  is  needed,  interpolation  must  be 
performed  on  the  range  compressed  signal  in  the  range  dimension.  The 
range  interpolation  can  be  done  in  the  azimuth  frequency  domain,  prior  to 
azimuth  feat  convolution  compression,  or  in  the  azimuth  spatial  domain, 
after  azimuth  compression.  In  the  former  approach,  the  range  curvature 
can  bo  fully  compensated,  whereas  in  the  latter  approach  only  a  piecewise 
compensation  can  be  done.  The  reason  for  this  is  as  for  lows:  once  the 
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azimuth  signal  is  compressed,  the  curvature  phase  is  transformed  into 
broadening  of  the  range  compressed  pulse,  and  the  trace  of  the  range  cur¬ 
vature  is  lost.  The  only  way  to  correct  for  curvature  after  the  azimuth 
signal  is  compressed  (i.e.  the  post-compression  scheme  previously  men¬ 
tioned),  is  to  use  sub-apertures,  each  of  which  produces  a  pulse  which  is 
minimally  distorted,  but  is  displaced  in  range  from  the  centre  of  the 
range  cell.  The  displacement  is  equal  to  the  amount  of  range  curvature  at 
the  sub-aperture  centre.  Range  shifting  of  the  sub- apertures,  by  means  of 
interpolation,  can  be  used  to  remove  the  curvature.  Such  interpolation  is 
equivalent  to  performing  a  piecewise  correction  along  a  single  azimuth 
array  in  the  frequency  domain.  This  is  shown  in  Figure  12.  This  piece- 
wise  curvature  compensation  is  incomplete  and  causes  paired-echoes  in  the 
final  image,  see  Sections  6  and  7  for  further  details. 


Fig.  12  -  Look  extraction  with  range  curvature.  (a)  Signal  history 
cf  a  point  target.  <b)  Look  extraction  with  piecewise 
range  curvature  correction  (only  the  positive  frequency 
axl s  is  shown ) . 


Despite  these  paired-echoes,  post-compression  interpolation  is 
preferred  in  the  present  scheme.  This  is  because  there  is  insufficient 
computer  memory  to  do  the  cross-row  interpolation  operations  on  the  long 
rows  needed  for  full  aperture  correction  in  the  range,  azimuth-frequency 
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domain.  In  this  case,  the  shorter  rows  required  for  the  piecewise 
correction  of  the  sub-apertures  make  the  post-compression,  range-azimuth 
spatial  co-ordinate  correction  technique  more  attractive. 

The  range  curvature  profile  ij>(t),  given  in  (43),  will  also  be 
approximated  in  a  piecewise  fashion,  namely. 


gq(ti,a)  *  *[tA  -  -  a2s2]  e 


"^-[■oq  +  »lq(3-Sq>  +  a2q<3"aq)2' 


2  2 

over  the  sub-aperture,  where  t*  *  —  a2s  near  the  look  centre.  >p(t)  now 
represents  a  coupled,  slanted qpro?ecti§n  of  the  compressed  range  profile 
onto  the  uncompressed  azimuth  profile.  This  coupling  only  affects  the 
spectrum  magnitude,  leaving  the  phase  unaffected.  Por  ease  of 
illustration,  <|>(t)  will  be  dropped,  i.e.,  it  will  be  assumed  the  curvature 
has  been  properly  corrected .  Thus,  (4 9)  becomes. 


9a(a)  *  « 


""}”■[  *oq  +  alq(3-3q)  +  a2q(s-Sq)2] 


The  phase  of  the  original  signal  is 


-[aQ  +  a^s  +  a2s2  +  a3S3] 


g(s)  »  e 


This  signal  is  ratified  by  removing  the  phase  term  -  — -  a^s,  during  data 
acquisition,  i.e.  the  recorded  phase  signal  is 


g(s)  «  g(s)  e 


+  a2s2  +  ajs3) 


It  now  remains  to  define  the  matched  filter  for  this  signal. 
Setting  g(s)agq(s)  over  the  sub-aperture  q,  we  get 


ga(s)  =>g(s)  ■  g(s)  e 


*  gq ( s )  e 


“^~~Laoq  +  alq(s~sq)  +  a2q(s_sq)2]  al3 


where  gQ(s)  is  a  quadratic  piecewise  approximation  to  g(s)  over  the 
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sub-aperture.  Equation  (5.17)  can 
"^[(&oq-alq3q)  + 

gq<s>  -  e 

and  the  matched  filter,  m(s),  can 
i.e. , 


be  regrouped  as, 

<alq-  aL)s  +  a2q(s-3q)2] 

,  (54) 

be  derived  from  this  regrouped 


gq(s) , 


ro( s )  ■  g  (-s) , 

( *oq  -  alqV  *  (alq"al>3  +  a2q( “3~3q ) 2 > 

■  e 

(55) 


The  Fourier  transform  of  m(s)  is 


MIC) 


P{m<a) } , 


1 ^aoq“alq3q5 
e 


ho  -*■  j2irsqh0 

e  4a2q 


(56) 


where  hf, 


(alq-ai). 


The  amplitude  of  M(f)  is  normalized  to 


unity,  and  the  small  amplitude  variations  from  look:  to  look  are  ignored. 
By  completing  the  square  of  the  phase  of  the  second  exponent,  (5.20)  can 
be  rewritten  as 


^•(*oq  -  alqaq  +  a2qaq>  — h2 

M(f)  -  e  A  e  4az<l  ,  (57) 

where  h  »  f  +  —  [(aiq  -  at)  -  2aqa2q]. 

A 

The  q  dependent  phase  terms  ensure  phase  continuity  at  the  boundaries 
between  adjacent  looks  in  the  concatenated  spectrum.  The  freouency 
dependent  phase  terras  cancel  the  quadratic  and  linear  phase  factors  in  the 
original  data.  The  above  matched  filter  can  be  applied  directly  onto  the 
Fourier  transformed  data  via  fast  convolution. 


6.  OOHPOTKR  SIMULATIONS  OF  OJHKHBff  SUB -APERTURE  PROCESSING  KITH 

PIECEWISE  CORRECTION 

In  this  section  the  digital  form  of  the  received  signal  from  a  unit 
amplitude  point  target,  and  the  form  of  the  continuously  straightened 
spectrum  of  the  signal  are  described;  the  effect  of  the  cubic  phase  term 
is  examined;  and  the  results  of  a  computer  simulation  of  piecewise 
correction,  i.e.  sub-aperture  processing,  are  presented. 

It  was  assumed  in  the  simulation  that  a  Hamming  window  was  applied 
to  the  range  frequency  data  prior  to  range  compression,  and  that  the  slope 
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of  range  curvature  versus  azimuth  position  in  the  aperture  was  small  at 
either  end  of  the  aperture.  It  was  also  assumed  that  the  linear  component 
of  range  migration  was  removed  during  the  data  acquisition  operation. 

6 . 1  Fon  of  the  Received  Signal . 

The  digital  form  of  the  received  signal  is 


,4it 


G(m,n)  -  R(m,n)  e  * 


(m) 


where  m,n  are  the  running  indices  for  azimuth  and  range 
R(m,n)  is  the  amplitude  profile  of  the  recorded  signal,  and 


r  (m) 


/ 


2  2 

(mASA)  +  rF  -  2mASArpcosnp  -  ai  mASA 


2  2 

(mASA)  +  rp  -  2mASArpCOsnp  +  mASAcosnp, 


respectively. 


(58) 


where  ASA  is  the  azimuth  sampling  interval  in  metres,  and 


m 


where  Ls  is  the  synthetic-aperture  length  in  metres. 


The  last  term  in  (58)  accounts  for  the  cancellation,  during  data 
acquisition,  of  the  linear  tern  mASAcosnp.  R(m,n)  is  the  oompressed 
range  profile,  which  runs  along  the  range  curvature.  It  is  given  by  the 
Fourier  transform  of  the  Hamming  window  (4],  i.e.. 


R(m,n)  ■  D i  +  0.426(D2  +  D3), 
where  »  sine', — ), 


(59) 


IT 

Nw 

TT 

**W 


X 

i) 

*  sinc(—  - 

1) 

Yr 

1 ) 

»  sinc(-—  + 

1) 

cr 

and  is  the  range  data  length  for  compression,  which  is  the  same  as 
the  window  length.  The  peak  amplitude  of  R(m,n)  has  been  normalized  to 
unity.  For  large  the  phase  terras,  in  D2  and  D3  are 

negligible.  The  variable  x  is  given  by 


x  ■  n  -  Nq 


r(m) 

w 


where  ASp  is  the  range  sampling  interval  in  metres,  is  the 


over-sampling  factor  in  range,  and  NQ  is  an  arbitrary  offset  for 
positioning  the  data  in  the  middle  of  the  two-dimensional  computer 
memory.  A  profile  of  a  compressed  range  pulse  is  shown  in  Figure  13.  To 
obtain  this  figure,  a  set  of  typical  operating  parameters  was  used  (see 
Table  1). 


OVEH-SAMPLING  FACTOR  «  1.1 
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MLg.  13  —  QaprMsad  range  poise  with 


weighting. 


6.2  CootimKnxa  Range  Curvature  Compensation 

As  previously  described,  one  way  to  remedy  the  energy  spill-over 
effect  caused  by  range  curvature  is  to  straighten  out  the  curvature  by 
continuous  correction  in  the  range-spatial,  azimuth- frequency  domain.  An 
azimuth  spectrum,  straightened  by  this  method  is  shown  in  Figure  14.  It 
is  virtually  identical  to  what  would  be  obtained  if  there  were  no 
curvature.  The  interpolator  used  in  the  straightening  process  was  a 
modified  four-point  sine  function  [5]  given  by 

W( x )  -  -3.ln(™)  (1  -  ill),  (60) 

7TX  16 

where  x  -  i  ♦  INTEGER Uj )-l\,  i  -  -1,  0,  1,  2,  (61) 


and  -  - — —  m 

8CTa2AsjiNT 


NT  Nrp 

t  n,  —  ■  x  m  x  • , 
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TABLE  1 

Processing  Parameters  for  Computer  Simulation  Examples 
Parameter 


Stant  Range  (rp)  100  km 

Wavelength  (X)  0.0321  m 

-3dB  azimuth  resolution  (p)  1.358  m 

Azimuth  sampling  interval  (ASA)  1.617  m 

Azimuth  over-sampling  factor  (yA)  1.18 

Range  sampling  intervals  (ASp)  0.75  m 

Range  sampling  period  (T)  5  nS 

Range  over-sampling  factor  (yR)  1.1 

Squint  angle  (np)  30* 

Range  curvature  at  the  end  of  the  full  2.33 

aperture 


N.B.  For  examples  with  squint  angles  (np)  equal  to  15*  and  6*,  A SA 
will  be  changed  accordingly,  all  other  parameters  are  kept 
constant. 


Pig.  14  -  Axisvth  spectrum  continuous  range  curvature 
ooopensat  ion . 

where  tip  is  total  number  of  data  points  used  in  processing.  Equation 
(62)  is  the  digital  form  of  (23)  with  f  ■  m/NTASA.  The  straightened 
spectrum  shown  in  Pigure  14  was  then  windowed  to  suppress  spectral  leakage 
and  multiplied  by  a  matched  filter.  The  windowed  spectrum  is  shown  in 
Pigure  15,  and  the  final  compressed  pul3e  is  3hovn  in  Figure  16. 
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SAMPLE  NUMBER 


rig.  16  -  (b)  sooa  in. 

6.3  Coble  Phase  Ten  Effect 

Aside  from  the  range  curvature  problem,  there  may  be  a  problem 
caused  by  the  cubic  phase  term.  If  the  cubic  phase  term  is  significant 
relative  to  the  quadratic  phase  term,  as  in  the  highly  squinted,  high 
resolution  case,  a  phase  correction  must  be  made  in  addition  to  the  range 
curvature  correction.  An  example  of  such  a  case  is  shown  in  Figure  17. 
In  this  situation  there  was  a  cubic  phase  error  of  1102*  at  the  end  of  the 
aperture.  This  cubic  phase  error  can  be  eliminated  by  adding  an 
appropriate  cubic  phase  term  in  the  matched  filter,  or  by  using  the 
sub-aperture  processing  technique. 

6 . 4  Sob-aperture  Processing  Computer  Simulations 

In  sub-aperture  processing,  the  spectrum  is  broken  up  into  small 
pieces  or  sub-apertures,  which  are  proportionally  shorter  in  length,  and 
can  have  their  FM  rates  fine-tuned  to  that  of  the  corresponding 
sub-aperture  data.  The  form  of  the  frequency  domain  matched  filter  is 

,  2,  2 

(aoq~a  lq3q+a2q3q^  ->—v 

A  4a2q, 


Min)  =  e 


e 


(63) 
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where,  v  ■  n  +  £NTA£.A  ((ajg-a^)  -  23^32^]; 

X 

Nt  is  the  FFT  length; 

Asa  is  the  azimuth  sampling  interval  in  metres; 

Asq  is  the  spatial  location  of  the  look  centre,  given  by  (55); 

f  .  1 

n  ■  — ;  and  Af  •  - . 

Af  nTASA 

Equation  (63)  is  the  digital  counterpart  of  (57). 
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Fig.  17  -  Compressed  pulse  with  cubic  phase  on compensated. 


Each  of  the  sub-apertures  is  filtered  with  its  own  M(n)  and  then 
the  spectra  for  the  sub-apertures  are  coherently  summed  by  concatenating 
them  together.  An  example  of  a  coherently  summed  spectrum  with  eight 
looks  is  shown  in  Figure  18  where  the  left  portion  of  the  spectrum  belongs 
to  negative  frequencies.  Prominent  gaps  are  evident  between  the  adjacent 
sub-aperture  spectra.  As  the  frequency  increases,  the  gaps  go  deeper. 
This  is  because  the  coupling  between  range  and  azimuth  increases  at  the 
ends  of  the  aperture,  i.e.,  where  the  range  curvature  slope  is  steepest. 

If  a  Hamming  window,  with  the  same  bandwidth  as  that  of  the 
original  spectrum  is  applied  to  the  concatenated  spectrum  shown  in  Figure 
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rig.  IS  -  ninth  spectra  with  piecewise  range  curvature 
condensation  (squint  angle  np  ”  30*). 
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Fig.  19  -  SpeetniM  of  Figure  18  with  Hanning  weighting. 
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19,  the  spectral  leakage  is  suppressed.  In  addition  and  perhaps  more 
importantly,  the  amplitude  of  the  gaps  in  the  envelope  are  equalized.  The 
gaps  modulate  the  envelope  intensity  and  cause  distortion  in  its  Pourisr 
transform.  The  equalization  effect  of  the  window  makes  the  envelope  modu¬ 
lation  almost  sinusoidal.  This  means  that  "paired-echoes"  will  be  pro¬ 
duced  in  the  other  domain  ( 1 ] .  Without  the  equalization  the  modulation  is 
non-sir.usoidal,  and  a  more  complex,  harder  to  analyze  distortion  is  pro¬ 
duced. 


An  output  pulse  from  one  sub-aperture  is  depicted  in  Figure  20. 
The  ma inlobe  width  is  wider  than  that  from  a  fully  compensated  full  aper¬ 
ture.  This  is  because  the  sub-aperture  bandwidth  is  proportionally 


SAMPLE  NUMBER 

Fig.  20  —  Compressed  poise  of  one  look. 


smaller  than  that  of  the  full  aperture.  Also,  there  are  prominent  side- 
lobes  .  This  is  to  be  expected  since  the  window  was  not  applied  to  the 
sub-aperture  spectrum.  Instead  the  window  will  be  applied  later  to  the 
concatenated  spectrum,  and  the  high  sidelobes  will  be  cancelled  out  at 
that  time . 

Thus  far  only  the  effects  of  amplitude  distortion  have  been  con¬ 
sidered.  Phase  errors  near  the  boundaries  cf  the  adjacent  looks  also 
cause  "paired-echo"  type  distortions.  The  phase  cancellation  during 
matched  filtering  is  incomplete  when  a  piecewise  correction  is  used.  This 
incomplete  phase  cancellation  causes  a  nearly  sinusoidal  residual  phase 
error  in  the  frequency  domain,  which  in  turn  causes  paired  echoes  in  the 
time  domain.  This  effect  is  serious  at  high  squint  angles  (e.g.  6*), 
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where  the  cubic  phase  term  (  —  833  )  becomes  very  large  and  cannot  be 

ignored.  A  plot  of  the  residual  phase  errors  for  our  present  example,  is 
shown  in  Figv  *e  21, 

In  order  to  obtain  the  final  compressed  pulse,  an  inverse  Fourier 
transform  of  the  windowed,  demodulated,  concatenated  spectrum  is  taken. 
The  result  is  shown  in  Figure  22 .  Except  for  the  presence  of  paired- 
echoes,  the  compressed  pulse  shape  is  identical  to  the  case  with  continu¬ 
ous  range  curvature  compensation.  As  has  already  been  mentioned,  these 
paired-echoes  are  the  products  of  envelope  and  residual  phase  distortions. 

Simulated  compressed  output  pulses  for  two  cases  with  high  scruint 
angles  (up  ■  15*,  6*)  are  shown  in  Figure  23  and  24.  As  can  be  seen  in 
the  figures,  the  paired-echo  magnitudes  increase  as  the  radar  becomes  more 
highly  squinted  (i.e.,  smaller  Up  is  used).  The  interrelationship  be¬ 
tween  various  radar  parameters  and  the  paired-echo  magnitudes  will  be  dis¬ 
cussed  in  the  next  section. 

7.  ERROR  ANALYSIS  OF  THE  EFFECTS  OP  PHASE  AND  AMPLITUDE  DISTORTIONS, 

CAUSED  BY  THE  PIET-ENISE  APPROXIMATION 

The  piecewise  correction  for  range  curvature  and  cubic  phase  intro¬ 
duces  phase  and  amplitude  errors  into  the  matched  filtering  process.  In 
fact,  exact  matching  occurs  only  at  the  look  centres.  Our  goal  here  is  to 
characterize,  with  respect  to  the  number  of  sub-apertures  and  the  other 
operating  parameters,  the  form  and  magnitude  of  each  of  the  distortions 
arising  from  these  abnormalities.  Because  an  exact  description  of  the 
errors  is  extremely  complicated  and  unnecessary,  a  combined  deterministic 
and  empirical  estimation  procedure  will  be  used  instead. 

In  this  section  the  form  of  the  magnitude  and  phase  errors  in  the 
frequency  domain,  will  be  examined.  Then,  models  will  be  developed  that 
allow  prediction  of  both  the  size  of  the  errors,  and  the  individual  and 
combined  paired  echo  levels  that  result  from  the  errors.  Next,  a  pro¬ 
cedure  for  estimating  the  Integrated  Sidelobe  Ratio  degradation,  caused  by 
the  errors,  will  be  described,  and  finally,  the  range  independence  of  a 
given  set  of  piecewise  corrections  will  be  examined. 

7 . 1  Spectrum  Envelope  Distortion 

In  the  following,  we  will  describe  a  simplified  approach  for  deter¬ 
mining  the  gap  size  between  adjacent  looks  in  the  frequency  domain,  as  a 
function  of  look  centre  location,  and  then  we  wi.1  1  relate  the  gap  sizes  to 
the  paired-echo  magnitudes  in  the  spatial  domain. 

When  post-compression  interpolation  is  used,  the  coupling  between 
range  and  azimuth  prevents  the  complete  correction  of  the  curvature,  and 
gaps  are  left  in  the  spectrum.  The  height  and  width  of  the  gap,  associ¬ 
ated  with  each  sub-aperture,  depend  on  the  compressed  pulse  width  in  range 
and  the  local  slope  of  the  range  curvature  versus  azimuth  position  curve. 
The  slope  of  the  range  curvature  function  in  the  frequency  domain  is  given 
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Fig.  22  -  Ccrap teased  pulse  with  piecewise  range  curvature 
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by  the  derivative  of  the  amount  of  range  curvature  (J^)  as  given  in  (62) 
with  respect  azimuth  space  (m) ,  i.e., 

slope  *  —  ” - - -  m,  (64) 

dm  4CTa2ASANT 

where  the  slope  is  expressed  in  number  of  range  cells  shifted  per  azimuth 
cell.  The  parameters  i  and  m  are  the  running  indices  in  range  and  azimuth 
respectively.  This  slope  is  used  to  project  the  compressed  range  profile 
onto  the  azimuth  spectrum.  It  is  assumed  that  a  Hamming  window  is  used  in 
range  compression,  and  that  therefore  the  compressed  range  profile  is  ap¬ 
proximately  the  shape  of  the  Fourier  transform  of  the  Hamming  window,  (see 
(59)). 


RANGE 


Fig.  25  -  Ranger- as inarth  coupling. 

The  coupled  range-azimuth  signal  for  one  sub-aperture  is  shown  in 
Figure  25.  It  is  assumed  that  the  aub-aperture  is  short,  that  the  quad¬ 
ratic  component  of  the  curvature  is  negligible,  and  that  the  slope 

(-^)  of  the  sub-aperture  segment  s  so  small,  and  so  constant  over  the 

3ub~aperture,  that  0(ra)  in  Figure  25  can  be  approximated  as: 

6(m)  s  tanO(m)  ■  —  <■  - - - - —  m*  (65) 

dm  4CTa  2AS^N<p 

With  this  approximation,  the  normalized  projection  of  the  compressed  range 
profile  onto  the  azimuth  frequency  axis  becomes, 

W(mq,p)  -  Ej  (Mq  ,p)  +  0.426[E2<mg,p)+E3(niq,p)  ] ,  (66) 

where, 

E;(mq,p)  -  ainc  (pk), 

E2(mg,p)  ■  sine  (pk-J), 

E3(mq,p)  “  sine  (pk+1). 
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k  -  — 

Yr  dm 
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m-nig  4CTa2AS^Nij'rR 
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_a 

Af 


p  ■  ro-Dig , 

and  fq  is  the  location  of  the  centre  of  look  q.  mq  denotes  the  centre 
of  the  qth  look  In  the  m  (global  frequency)  dimension,  p  is  a  new  local 
frequency  dimension,  whose  origin  is  at  the  centre  of  the  look. 
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Pig.  26  •*  Look  extraction  using  a  low  pass  filter. 


As  shown  in  Figure  26  the  desired  sub-aperture  is  extracted,  by 
means  of  an  FIR  filter,  from  an  appropriately  rotated  version  of  the 
baseband  spectrum.  If  the  FIR  filter  shape  is  approximately  rectangular, 

Nr 

the  spectral  height  at  the  filter  edge  (p  -  — )  is 


hq  -  W(mq,p)  -  Ex(mq,  +  0.426  [E2(Mq,  -~) 

Nr.  (67) 

+  E3(raq,  )], 

where  Np  is  the  number  of  frequency  cells  covered  by  the  bandwidth  of 
the  extraction  filter. 

For  the  looks  extracted  from  the  ends  of  the  spectrum,  the  slope  of 
the  range  curvature  .<  s  steeper,  and  the  area  of  the  projected  range 
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profile  i.a  narrower.  As  a  result,  the  gap  between  adjacent  sub-apertures 
is  widened  and  deepened.  In  other  words,  the  envelope  distortion  worsens 
towards  either  end  of  the  full  aperture.  A  typical  coherently  summed 
spectrum  is  sketched  in  Figure  27.  As  shown  in  the  figure,  the  minima  can 
be  connected  by  a  dotted  line  described  by, 

h(m)  -  Fj(m)  +  0.426  [F2{ni)  +  ^3 (M) ] ,  (68) 


where. 


and. 


F i (m)  -  sine  (Ko) , 

P2 (m)  ■  sine  (Rm-1 ) , 

P2 (m)  -  sine  (Km+1 ) , 

K  - - —  . 

8CTa2AS^YR 


(m> 

Fig.  27  -  A  typical  coherent  sunaed  spectrum.  The  dotted  line 
connecting  all  the  Minima  is  described  by  Bq«,  68. 


As  mentioned  previously,  Hamming  weighting  is  normally  applied  to  the 
concatenated  spectrum  to  suppress  the  spectral  Leakage  that  occurs  during 
the  calculation  of  the  inverse  Fourier  transform.  The  weighting  function 
has  the  name  bandwidth  as  the  summed  spectrum.  The  weighted  spectrum  is 
sketched  in  Figure  28 . 

Besides  suppressing  spectral  leakage,  the  weighting  also  helps 
reduce  the  size  of  the  larger  gaps  at  the  ends  of  the  concatenated 
spectrum.  The  gaps  are  now  more  or  less  equalized  over  the  entire 
aperture.  As  shown  in  Figure  2 V,  the  weighted  gaps  can  be  thought  of  as 
ripples  superimposed  on  an  ideal  smooth  spectrum.  The  amplitude  of  the 
ripples  at  a  distance  of  a  quarter  of  the  bandwidth  from  the  aperture 
centre  is  estimated  to  be. 
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Fig.  28  -  Spectral  energy  distribution. 
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Fig.  29  -  Save lope  distortion  (assune  sinusoidal  ripples)  for 
coherent  snlti-look  spectra*. 


where  h(m)  is  given  by  (29)  and  the  0.54  numerical  factor  accounts  for  the 
Hamming  weighting.  This  ripple  amplitude  is  taken  to  be  representative  of 
the  ripple  across  the  whole  aperture. 
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The  ripples  shown  in  Figure  28  exhibit  a  periodic  sinusoidal-like 
waveform.  Assuming  for  now  that  the  ripples  are  purely  sinusoidal  (which 
is  not  really  true),  the  inverse  Fourier  transform,  as  shown  in  Figure  29, 
will  give  paired-echoes  of  magnitude  hj/2,  on  either  side  of  the  mainlobe 
of  height  approximately  equal  to  0.54. 


(70) 

ftp 

where  m’  -  -  . 

*tr 


Onfortunately,  the  computer  simulations  did  not  produce 
paired-echoes  of  the  predicted  magnitude.  One  explanation  for  thi3 
disagreement  is  that  the  ripples  on  the  sunned  spectrum  are  not  purely 
sinusoidal.  If  this  is  the  case,  the  first  order  harmonic  coefficient 
should  fairly  cloaaly  represent  the  most  dominant  paired-echo.  In  the 
absence  of  an  analytical  expression  for  the  ripple,  it  is  impossible  to 
prove  this  implication.  Instead,  a  correction  factor  was  established, 
which,  when  used  in  the  simulation,  caused  paired-echoes  of  the  correct 
magnitude  to  be  produced.  It  was  found  that  if  a  correction  factor  of 
one-half  was  used  In  (70),  the  results  were  reasonable  in  most  cases. 
With  a  correction  factor  of  ona-half,  (70)  becomes 


main  peak 


f 1-h' ,  for  more  than  two  looks; 

I  8 

1-h ' ,  for  two  looks. 

12 

(71  ) 


A  special  provision  is  made  for  the  two-look  case,  i.e.,  the  correction 
factor  is  set  to  1/3  instead  of  1/2.  This  is  required  because  the 
envelope  ripple  deviates  more  strongly  from  a  sinusoidal  shape  in  the 
two- look  case  than  it  does  in  the  other  cases. 


These  results  are  consistent  with  the  fact  that  as  the  ripple  shape 
deviates  more  and  more  from  the  sinusoidal,  more  energy  is  displaced  into 
the  higher  order  harmonics,  and  thus  the  magnitude  of  the  dominant 
harmonic  decreases. 
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With  changes  in  the  operating  parameters,  the  general  shape  of  the 
ripples  is  preserved,  and  the  relative  values  of  Fourier  coefficients  do 
not  change.  Equation  (71)  provides  a  good  estimate  of  the  paired-echo 
magnitudes  when  the  magnitudes  are  between  -20dB  to  -35dB.  For  ratios 
greater  than  -20dB,  the  aforementioned  linear  approximation  assumptions 
break  down  gradually,  with  the  estimated  ratio  increasing  more  rapidly 
than  the  actual  ratio.  For  ratios  less  than  -35dB,  (71 )  also  breaks  down, 
this  time  because  the  pair-echo  magnitude  becomes  so  low  that  it  merges 
with  the  shoulder  of  the  main  peak. 


The  location  of  the  dominant  paired-echo,  which  is  related  to  the 
number  of  looks  involved,  is  given  by 


nDe  =  —  cells, 
P 


(72) 


where  is  the  number  of  looks  and  p  is  the  ratio  of  the  overall 

extracted  spectrum  to  the  Nyguist  bandwidth.  In  most  cases. 


0.9, 


(73) 


and  thus 

npe 


cells. 


(74) 


In  the  above  analysis  and  simulations,  envelope  distortion  was 
conaidered  to  be  the  sola  contributor  to  paired-echoes,  and  the  phase  was 
assumed  to  be  perfectly  matched.  In  the  next  section,  the  effects  of 
phase  errors  are  examined. 


7.2  Residual  Phase  Error 


The  residual  phase  error  originates  from  a  phase  mismatch  between 
the  piecewise  quadratic  matched  filter  and  the  data.  This  error  causes 
paired- echoes . 

The  phase  error,  between  the  recorded  data  and  the  matched  filter, 
at  s“0 ,  is  taken  to  be  representative  of  the  phase  error  ripple  across  the 
entire  aperture.  This  representative  phase  error  is  obtained  by 
subtracting  the  phase  in  (4  5)  from  (47),  and  then  setting  s«0,  and  g**1 , 
i  .e. 


Ae  - 


4* 

r 


«ol  ~  *1131 


+ 


*21°  1 


2 


*o 


(75) 


Baaed  on  the  correspondence  between  the  spatial  and  frequency  domains 
described  by  (21)  it  is  possible  to  find  an  fi  for  every  s^,  such  that  the 
phase  iii  the  freauency  domain  at  fj  equals  the  phase  in  the  spatial  domain 
at  s  i .  If  the  principle  of  stationary  phase  holds,  it  also  follows  that 
the  maximum  phase  error  in  the  frequency  domain  is  equal  to  the  maximum 
phase  error  in  the  spatial  domain.  Therefore  the  residual  phase  error. 
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A0,  given  by  (75),  is  equal  to  the  phase  error  at  the  boundary  of  the 
first  look  in  the  frequency  domain. 

It  was  assumed,  and  confirmed  later  by  computer  simulations,  that 
the  peak  residual  phase  errors  at  the  end  of  the  sub-apertures  are  approx¬ 
imately  equal. 

A  sketch  of  the  typical  residual  phase  error  in  a  coherently  summed 
spectrum  is  shown  in  Figure  30.  The  residual  phase  error  is  periodic  in 
pattern.  Because  of  the  overlap  (18  to  21%)  of  the  transition  bands  of 
the  look  extraction  filters,  the  residual  phase  errors  at  the  look  bound¬ 
aries  cancel  each  other  out.  This  cancellation  reduces  the  peak  error 
(A0O) ,  to  about  2/3  of  A9,  i.a. 

A9p  *  1  A0 

“  ~  |  *ol  *  *llal  +  a213l2  ~  |-  <76) 


If  the  phase  error  is  purely  sinusoidal,  with  amplitude  A0p,  where  A0p 
<  0.4  radian,  then  the  paired-echo  magnitude  is  given  by  a  Bessel  function 
approximation  [ 1 ] ,  such  that 


paired-echo  due  to  phase  error 
main  peak 


(77) 


RESIDUAL  PHASE  ERROR 
(ASSUME  SINUSOIDAL) 
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As  was  the  case  with  the  envelope  distortion,  the  phase  error  is 
not  exactly  sinusoidal  and  no  simple  solution  can  he  formulated.  However, 
the  size  of  the  dominant  first  harmonic  coeffit  ent  can  be  obtained  by 
computer  simulation.  When  this  is  done,  and  20%  transition  band  is 

assumed,  the  magnitude  of  the  first  harmoni  j  out  to  be  only  about 

half  of  that  given  by  (77),  i.e., 

paired-echo  due  to  phase  error  __  A8P 
^  main  peak  4  (78) 

Since  the  residual  phase  error  in  the  frequency  domain  has  the  same  period 
as  that  of  the  envelope  distortion,  the  locations  of  the  paired-echoes, 
caused  by  the  phase  distortion,  coincide  with  the  locations  of  the  echoes 
produced  by  the  envelope  distortion. 

7.3  Overall  Paired-echo  Magnitude 

The  first  order  paired-echo  magnitudes  are  given  ^v  the  sum  and 
difference  [1]  of  that  contributed  by  envelope  distortion 

(70)  and  residual  phase  error  (78),  i.e., 

overall  paired-echo  =■  IPmJtPpj*  (79) 

main  peak 

7.4  Degradation  of  Integrated  Si  delobe  Ratio  doe  to  Paired-echoes 

For  a  heavily  weighted  spectrum  (e.g.  Hamming),  in  the  absence  of 
additional  phase  errors,  e.g.  FM  rate  mismatch  or  motion  errors,  the 

integrated  sidelobe  level  (everything  outside  the  first  minima  from  the 
main  peak)  is  composed  mainly  of  energy  contributed  by  the  paired-echoes, 
|PmiPp|*  integrated  sidelobe  ratio  (ISLR)  can  be  approximated 

as 


XSLR  *  10  logi0[  (pnj  +  pp)2  t  (Pm-Pp)2] 

(80) 

m  10  +  Pp  )]. 

7  .  S  Range  Independence  of  the  Piecewise  Solution 

The  resultant  piecewise  solution  can  be  regarded  as  a  one- look 
matched  filtering  process  with  a  distorted  frequency  transfer  function. 
The  amplitude  versus  frequency  response  is  modulated  with  sinusiodal-like 
ripples  in  its  passband,  while  the  phase  versus  frequency  response  is 
composed  of  ripple- like  residual  phase  errors  superimposed  on  the  usual 
linear  and  quadratic  phase  cont  ributions.  In  both  cases,  the  distortions 
are  local  in  nature  since  they  arise  from  mismatches  within  a  look. 

In  a  normal  one-look  azimuth  matched  filtering  process,  the  filter 
is  matched  to  the  azimuth  signal  returned  from  a  radar  point  target 
located  at  the  centre  of  the  processing  swath.  This  matched  filter  is 
then  applied  to  all  the  data  in  that  processing  swath.  The  quadratic  or 
cubic  phase  mismatch  between  the  data  and  the  filter  increases  towards  the 
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edges  of  the  swath.  The  width  of  a  swath  is  determined  by  the  maximum 
amount  of  phase  mismatch  tolerable,  at  the  end  of  the  processing 
aperture.  This  swath  width  is  known  as  the  depth  of  focus  (2], 

During  the  application  of  the  azimuth  nulti-look,  frequency  domain, 
matched  filter  to  data  at  different  ranges  across  the  swath,  the  compress¬ 
ed  frequency  domain  signal  will  suffer  both  local  and  global  phase  mis¬ 
matches,  as  well  as  envelope  distortion.  The  overall  distortion  in  the 
frequency  domain  signal  will  be 

j*e(f)  +  j4>g(f) 

e(f)  -  A(f)  e 

%<f) 

■  [a( f )  e  ]  e 

In  the  spatial  domain  the  error  will  be 

e(s)  -  P“ 1 { e ( f ) }  -F-1{A(f)e  }*F~T (e  }, 

(8',) 


where, 

A C f )  is  the  envelope  ripple  as  a  function  of  frequency; 

<Jie(f)  is  the  local  residual  phate  error  as  a  function  of 
frequency;  and 

4>g(f)  is  the  global  phase  mismatch  an  a  function  of 
frequency. 

In  (81),  the  local  phase  error  can  be  combined  with  the  envelope 
distortion  and  inverse  Fourier  transformed.  The  global  phase  can  be  in¬ 
verse  Fourier  transformed  separately.  The  first  inverse  Fourier  transform 
gives  rise  to  paired-echoes,  whereas  the  second  degrades  the  integrated 
aidelobe  ratio.  The  convolution  of  the  two  inverse  Fourier  transforms 
describes  the  total  effect.  Thus,  in  evaluating  the  paired-echo  magni¬ 
tudes  only  the  localdistortions  need  to  be  considered,  whereas  in  evalu¬ 
ating  the  degradation  in  the  integrated  aidelobe  ratio  (separate  from  that 
contributed  by  the  paired-echoes),  only  the  global  phase  mismatch  needs  to 
be  considered. 

8  .  DERAMPING  TBCHH1QPBS 

This  section  is  an  overview  of  alternate  farms  of  pulse  compression 
foe  a  linear  FM  signal.  Instead  of  matched  filtering,  or  fast  convolu¬ 
tion,  a  deramping  (dechirping),  spectral  analysis  approach  can  be  used. 
In  the  conventional  deramping  method,  the  quadratic  phase  term  in  the 
incoming  radar  return  signal  is  cancelled  by  multiplication  with  a 
reference  signal,  which  has  a  phase  term  that  is  the  conjugal*  of  the 
returned  signal.  The  process  is  shown  diagrainatically  in  Figure  The 

deramped  signal  becomes  a  CW  signal,  which  has  a  frequency,  flf  directly 
proportional  to  the  time  separation,  A,  between  the  received  signal  and 
the  reference.  The  CW  signal  is  then  Fourier  transformed  to  give  a 
compressed  pulsa  as  shown  in  Figure  31(c).  The  location  of  the  pulse  is 
directly  proportional  to  the  deramped  CW  frequency,  f^.  As  a  result,  the 
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target  is  resolved  at  its  location  of  arrival,  A.  Figure  32  shows  a 
multi-target  environment,  in  which  three  targets  are  resolved  at  their 
respective  locations. 

The  width  and  magnitude  of  the  compressed  pulse  depend  on  the 
duration  of  the  overlap  (Te)  between  the  reference  and  the  received 
signal,  as  shown  in  Figures  32(a)  and  (b).  For  targets  which  lie  further 
away  from  the  reference  (i.e.  smaller  Te),  the  derawped  signals  suffer 
more  energy  loss,  and  therefore,  their  pulse  width  and  magnitude  are 
degraded,  a a  shown  in  Figure  32(c).  In  practice,  there  exists  a  threshold 
above  which  the  amount  of  energy  loss  is  intolerable.  In  order  to  stay 
within  the  threshold,  it  is  necessary  to  overlap  reference  signals  so  as 
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Fig,  31  -  Conventional  deraaping  technique  for  a  point  target. 


to  decrease  the  time  separation  between  the  received  signal  and  the 
nearest  reference.  Such  a  scheme  reduces  the  amount  of  energy  loss,  but 
at.  the  expense  of  having  the  data  repeatedly  deramped  by  nearby  reference 
signals.  This  is  illustrated  in  Figure  33.  The  repetitions  produce 
redundant  data  with  energy  loss  above  the  tolerable  threshold  Only  data 
in  the  vicinity  of  the  centre  of  a  reference  ramp  are  kept.  The  redundant 
data  are  discarded.  The  valid  data  derived  from  successive  references  are 
concatenated  to  form  the  final  image. 


Pig.  32  -  Conventional  derating  technique-  far  milt i- targets. 


pig.  33  -  Overlap  reference  ra*p*  to  redace  collapsing  loss  . 


In  unis  section,  we  will  overview  two  pulse  compression  techniques 
based  on  the  aforementioned  deramping  approach,  nt>-v_iy,  the  Spectral 
Analysis  [ 1 0 ]  and  the  Step  Transform  (6-9]  Methods.  Compared  to  the  basic 
de-ramping  approach  both  methods  allow  a  reduction  of  tne  energy  loss 
and/or  an  improvement  in  the.-  processing  efficiency.  The  following 
description  fo-ms  a  fremewor-.  tor  the  discussion,  in  Section  9,  oir  the 
processing  efficiency  oi  *. he  sub-aperture  compression  techniaues. 
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8.1  Spectral  Analysis  Method  (9,10] 

The  spectral  analysis  method  attempts  to  reduce  the  energy  lost, 
due  to  aliasing,  during  the  deramping  process.  This  is  achieved  by 
utilizing  a  continuous,  non-overlapping  deramp  sawtooth  reference  with  an 
FM  rate  of  the  same  magnitude,  but  of  opposite  slope  and  longer  duration, 
than  the  signal  from  a  single  point  target.  In  addition,  the  phase  at  the 
end  of  one  sawtooth  is  it  and  that  at  the  beginning  of  the  next  sawtooth  is 

-it,  or  vice  versa  (i.e.  the  phase  at  the  boundary  is  circularly  continuous 

and  anti-symmetric).  A  signal  ramp  crossing  a  sawtooth  boundary  does  not 
change  its  deramped  CW  frequency  because  aliasing  causes  the  disjoint 
frequency  to  fold  back  in  the  spectrum  and  align  with  the  unaliased 

portion.  The  phase  is  continuous  at  the  boundary.  This  is  illustrated  m 

the  frequency  versus  time  diagram  shown  m  Figure  34(a)-(c).  It  i3 
therefore  possible  to  employ  a  Fourier  transform  to  extract  uhe  spectrum 


f 


Fig.  34  -  Ifee  spectral  analysis  method. 
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representing  the  point  target. 

For  multiple  targets,  superimposed  target  returns  can  be  uniquely 
resolved  by  the  deramping  -  Fourier  transform  process.  Figure  35(a)  shows 
a  parade  of  equally  spaced  target  returns  in  the  frequency  versus  time 
plane.  Figure  35(b)  shows  the  corresponding  deramping  reference  signal. 
Note  that  it  has  opposite  FM  slope  as  the  target  signals.  Upon  mixing  the 
target  signals  with  the  reference,  the  individual  FM  signals  are  trans¬ 
lated  into  CW  signals,  with  their  frequencies  corresponding  to  their 
positions  relative  to  the  reference.  The  resultant  deramped  signal 
ensemble  13  depicted  in  Figure  35(c). 


Fig.  35  -  Deramping  process  via  the  spectral  analysis  method. 


Signal  ramps  falling  within  a  period  of  the  sawtooth  reference  are 
grouped  into  a  parallelogram  in  the  frequency  versus  time  plane.  The 
oblique  boundaries  between  parallelograms  prohibit  Fourier  transformation 
of  the  full  aperture  data  set.  Processing  can  only  be  performed  with 
shorter  Fourier  transforms  that  prevent  data  from  the  adjacent  parallelo¬ 
gram  from  being  drawn  into  the  transform  and  getting  mixed  up  with  the 
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returns  that  belong  there.  It  should  be  noted  that  the  data  are  only 
directly  accessible  in  the  cime-domain  (i.e.  column-wise  in  Figure  35(c)), 
because  all  the  targets  are  superimposed.  The  frequency  dimension 
(vertical  axis)  is  accessible  through  either  the  DFT  or  filtering 
separations  on  the  time  domain  data.  As  already  mentioned,  one  way  to 
avoid  mixing  up  the  data  is  to  partition  the  parallelogram  into  narrow 
columns.  Each  column  is  centred  in  the  neighbourhood  where  the  target 
data  are  to  be  resolved.  Each  column  of  data  (actually  a  superimposed 
one-dimensional  time-domain  data  ensemble)  is  Fourier  transformed  along 
the  time  axis  to  resolve  the  target  ensemble.  A  data  column  will 
tnevitably  cross  a  parallelogram  boundary,  thus  mixing  up  data  between  the 
parallelograms.  Compressed  data  in  the  vicinity  of  a  boundary  are  invalid 
and  oust  be  discarded.  The  placement  of  one  processing  region  is  shown  in 
Figure  36. 
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Fig.  36  -  Placement  of  one  processing  region. 


Because  of  the  antenna  aperture  modulation,  target  signals  com¬ 
pressed  using  thir,  techriicjue  are  subject  to  amplitude  modulation.  The 
amplitude  nodulat.i  v  can  be  controlled  by  keeping  the  data  of  interest 
away  from  Dio  edge.  >.  ■  the  parallelogram. 


Targets  in  the  remaining  portions  of  the  parallelogram  can  be 
similarly  compressed  by  partitioning  another  processing  column  at  the 
other  end  of  the  parallelogram  as  shown  in  Pigure  37.  Compressed  data 
from  the  two  processing  regions  within  the  parallelogram  are  concatenated 
to  form  l.  continuous  image.  In  practice,  a  guard  band  of  fractional  width 
6  is  included  to  keep  the  valid  processed  data  away  from  the  oblique 
border.  As  mentioned  previously,  this  is  done  to  control  amplitude 
modulation.  A  typical  value  of  6  should  be  0.3  to  0.5  depending  on  the 
antenna  pattern.  The  output  of  an  N-point  DPT  computed  along  a  processing 
column  spans  the  entire  frequency  space  of  B  Hz.  When  weighting  is 
applied  to  the  N-point  DPT,  the  effective  DPT  length  on  (where  a  <  i)  is 
chosen  to  be  less  than  or  equal  to  the  -3dB  width  of  the  weighting 
function.  This  narrower  processing  column  is  then  used  to  calculate  the 
number  of  valid  output  points.  With  the  aid  of  Pigure  36,  the  number  of 
valid  points  is  found  to  be 

G  -  (l-0  -  £^)n  samples,  (02) 

B2 

where  k  is  the  azimuth  linear  PM  rate? 

B  is  the  reference  sawtooth  bandwidth;  and 

N  is  the  DPT  length,  which  is  fixed  for  all  the  processing 
columns. 


H2  samples 


Pig.  37  -  Placement  of  naccesaive  prccesing  regions  to  ensure 
image  continuity , 
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The  invalid  data,  discarded  because  oE  the  boundary  crossing  problem,  can 
be  regained  by  putting  another  processing  column  near  the  opposite  end  of 
the  parallelogram,  aa  shown  in  Figure  37.  This  second  processing  column 
has  the  same  width  (N  points)  as  the  first  one.  Its  position  is  deter¬ 
mined  by  the  frequency  at  which  the  first  column  is  cutoff.  The  second 
processing  column  encompasses  both  data  left  over  from  the  first 
processing  column  and  data  from  the  succeeding  parallelogram.  Part  of  the 
output  is  invalid  data  corresponding  to  the  frequencies  which  straddle  the 
oblique  boundary.  In  order  to  ensure  image  continuity,  the  second 
processing  column  is  located  such  that  the  beginning  of  the  region  of 
valid  output  samples  in  the  second  DFT  is  at  the  same  frequency  as  the  end 
of  the  corresponding  region  of  the  first  DFT.  This  is  shown  as  a 
horizontal  dotted  line  in  Figure  37.  Prom  simple  geometry,  t  e  required 
continuous  frequency  coverage  is  achieved  if  the  second  DFT  i.  begun  at 

S  -ia2-  (1-8  -  25*)  *1  -  (,-e)  -  »»  (83) 

N  k  B2  k  k 

input  samples  after  the  start  of  the  first  DFT. 

Since  N  is  not  equal  to  tha  width  of  a  parallelogram,  the 
processing  columns  are  not  synchronized  with  the  parallelogram 
boundaries.  As  the  processing  progresses,  the  boudnary  cuts  the 

processing  columns  at  different  frequencies,  therefore  the  locations  of 
the  valid  data  regions  must  be  updated  for  each  processing  column. 

The  above  analysis  assumes  that  two  processing  columns  are 

sufficient  to  resolve  the  entire  bardvidth  with  the  prescribed  resolution, 
but  that  is  not  necesarily  true.  Situations  may  occur  where  the  required 
processing  column  is  so  wide  and  0  (guard  band)  is  so  large  that  more  than 

two  processing  columns  are  required  within  a  single  parallelogram.  An 

example  of  such  a  situation  ia  shown  in  Figure  38 .  In  the  case  depicted 
in  the  figure,  none  of  the  processing  columns  produce  any  valid  data  in 
the  adjacent  parallelogram.  In  general,  this  happens  when 

G  <  —  ,  i.e. , 

2 


1-0 


or 


(34) 


If  a  substitution  involving  the  resolution  pa,  is  given  by 

1 . 4VB 

P  a  -  f 

oNk  (85) 

where  V  is  the  velocity  of  the  radar,  is  made  in  (84),  the  inequality  can 
be  written  as 


r 
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B  >  -  -  — —  .  (86) 

2  Bpa 

Note  in  (85)  that  N  is  restricted  to  integer  powers-of-two.  This  places 
limits  an  the  permissible  values  pa  may  assume. 
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38  -  PlaoMwnt  of  processing  regions  for  a  case 
with  large  guard  band. 


Major  drawbacks  of  this  method  include:  cubic  phase  errors  and 
range  curvature  cannot  be  compensated  easily,  it  is  difficult  to  achieve 
full  resolution  efficiently,  and  there  is  a  variation  of  signal-to-noise 
ratio  over  the  image  due  to  the  fact  that  different  sections  of  the  scene 
are  illuminated  by  different  portions  of  the  antenna  pattern. 

8.2  Step  Transform 

The  aim  of  the  step  transform  technique  is  to  reduce  the  energy 
loss  encountered  with  high  resolution  deramp  processing,  while  maintaining 
a  high  processing  efficiency  and  efficient  memory  utilization.  This  is 
done  by  conforming  the  processing  path  to  the  trend  of  the  signal  ramp, 
thereby  avoiding  the  problem  of  having  to  use  multiple  full  length,  highly 
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overlapped  reference  signals.  This  method  is  efficient  for  compressing 
linear  FM  signals  in  a  minimal  memory  space,  bat  it  lacks  on  inherent  cap¬ 
ability  for  accommodating  cubic  or  other  non-linear  phase  characteristics 
of  the  frequency  chirp.  In  the  next  section,  the  effect  of  the  cubic 
phase  term  will  be  investigated,  and  a  means  for  avoiding  image  degrada¬ 
tions  caused  by  cubic  errors,  suggested. 

In  che  step  transform,  the  multiple  references  in  the  frequency 
versus  time  plane  are  replaced  by  a  continuous  sawtooth  reference  func¬ 
tion,  as  shown  in  Figure  39.  The  sawtooth  reference  will  overlap  an  in¬ 
coming  signal  more  evenly,  regardless  of  its  time  of  arrival,  and  thus 


Pig.  39  -  Step  transform  processing  using  a  continuous 
sawtooth  reference  function. 


reduce  the  energy  loss  to  almost  zero.  Mixing  an  incoming  rami  with  the 
sawtooth  reference  gives  a  step-like  waveform,  in  the  frequency  versus 
time  plane,  as  shown  in  Figure  40.  The  horizontal  part  of  each  step  is  a 
CW  waveform.  The  vertical  parts  of  the  steps  are  of  equal  size  and  are 
spaced  evenly  in  time.  This  step-wise  partitioning  of  the  full  aperture 
return  signal  is  similar  to  the  sub-aperture  processing  described  in  the 
previous  sections.  In  the  matched  filtering  case,  the  look  extraction  was 
performed  by  a  combination  of  FIR  and  matched  filter  in  the  frequency 
dceutin.  Hence  each  extracted  look  consists  of  scenes  illuminated  by  the 
same  portion  of  the  antenna  pattern.  In  the  step  transform,  each  look 
extraction  is  performed  by  a  Hamming  weighted  conjugate  ramp  in  the  time 
doeoifi,  Hence  each  extracted  look  consists  of  scenes  illuminated  by  a 
different  part  of  the  antenna  pattern. 
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rig.  46  -  signal  alter  damped  By  a  oontlnooae  sawtooth 
reference  function. 


In  the  step  transform,  as  in  the  basic  deramping  techni<rue 
previous!,  ^c-L.ibed,  dl  ffeLtu.t  targets  will  generate  different  Doppler 
frequencies  depending  on  their  locations  relative  to  the  aub-ranrps.  If 
part  of  the  signal  data  for  several  different  targets  falls  into  the  same 
sub-ramp,  the  extracted  look  will  be  composed  of  an  ensemble  of  CW  signals 
whose  frequencies  are  directly  proportional  to  the  locations  of  the 
targets  on  the  sub-ramps.  The  frequency  produced  by  mixing  the  signal 
with  the  reference  will  also  be  directly  related  to  position  in  the 
azimuth  antenna  pattern.  Thus,  it  can  be  seen  that  different  targets  in 
the  aub-aperture  are  acquired  from  different  portions  of  the  antenna 
pattern.  This  is  the  basic  difference  between  the  matched  filter  and  the 
step  transform  approaches.  The  scenario  for  the  step  transform  is 
illustrated  in  Figure  41 .  The  image  formed  by  the  first  stage  of  the  step 
transform  is  of  low  resolution  because  only  a  small  fraction  of  the  target 
energy  is  contained  in  each  sub-aperture.  The  remaining  portions  of  the 
signal  energy  are  embedded  in  the  other  sub-apertures  that  together  make 
up  the  full  aperture.  In  order  to  regain  the  original  resolution,  the 
sub-apertures  must  be  added  coherently. 

The  case  of  an  input  signal  consisting  of  one  point  target  is  shown 
in  Figure  42 .  The  spectrum  associated  with  the  target  moves  to  a  higher 
frequency  as  it  climbs  each  step  up  the  frequency  staircase.  In  addition, 
there  is  a  step  increment  in  linear  phase,  which  depends  on  the  separation 
in  time  of  the  signal  and  the  reference  tooth.  These  two  properties  can 
be  used  to  regain  the  original  high  resolution  of  the  input  signal. 


59 


ONLY  THIS  COLUMN 
IS  SHOWN 


Pig.  41  -  Derawping  of  aoltipls  targets. 

The  frequency  step  in  a  sub-aperture  (see  Figure  41(b))  can  be 
segregated  from  the  superimposed  ensemble  data  by  utilizing  a  PFT  as  a 
filter  bank.  The  frequency  stepping  behaviour  necessitates  a  two- 
dimensional  processing  scheme  indexed  by  spectral  frequency  and  sub- 
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aperture  number,  aw  illustrated  in  Figure  42.  All  the  spectral  energy 
associated  with  a  single  target  can  be  recaptured  by  coherently  summing 
the  sub-apertures  along  a  diagonal  line.  Each  line  of  the  diagonally 
arranged  spectra  associated  with  a  point  target  contains  a  linear  phase 
whose  frequency  is  linearly  dependent  on  its  distance  from  the  nearest 
sub-aperture  on  the  time  axis.  The  form  of  this  phase  term  will  be 
derived  in  the  next  section.  Because  of  this  intimate  relationship 
between  local  frequency  and  time,  a  diagonal  Df*T  can  be  used  to  3ort  out 
the  fine  details  which  are  hidden  in  the  first  set  of  DFTs. 


r 


DIAGONAL 

DFT 

FOR  FINE 
RESOLUTION 


Fig.  42  -  Processing  matrix  after  the  first  set  of  DFTs. 


In  performing  the  first  DFT  on  each  sub-aperture,  windowing  must  be 
applied  in  order  to  suppress  sidelobes  and  spectral  leakage.  It  is  essen¬ 
tial  to  suppress  these  effects  so  that  energy  spill-over  onto  the  neigh¬ 
bouring  diagonals  is  minimized.  One  major  negative  aspect  of  the  window¬ 
ing  is  that  the  spectrum  is  broadened.  Since  the  first  DFT  is  regarded  as 
a  filter  bank,  a  broadened  spectrum  due  to  windowing  means  a  wider  band¬ 
pass  filter  element  (wider  them  one  cell).  The  widened  bandwidth  will 
undoubtedly  let  in  energy  from  the  adjacent  diagonals.  These  intruders, 
which  are  further  away  from  the  diagonal  under  consideration,  contain 
higher  frequency  components.  These  higher  frequency  components  exceed  the 
sampling  rate  along  the  diagonal.  If  a  DFT  is  taken  along  the  diagonal 
data  array,  the  intruding  components  will  fold  into  the  desired  region. 
In  order  to  avoid  aliasing,  the  sampling  frequency  along  the  diagonal 
direction  must  be  increased.  In  order  to  do  this,  the  dimension  of  the 


processing  matrix,  i.e.,  the  spacing  between  the  first  DFTs  must  be 
decreased.  Since  a  minimum  size  DFT  is  required  for  adequate  resolution 
from  the  first  processing  stage,  the  first  DFTs  become  overlapped  when  the 
spacing  is  decreased.  The  above  explanation  is  summarized  in  Figure  43 
[7].  A  block  diagram  of  the  step  transform  processing  is  shown  in  Figure 
44. 


Pig.  43  -  Tine  weighting  and  overlap  to  reduce  time  aide  lobes  and 
aliasing  [71. 

8.2.1  Formulation  of  the  Step  Transform 

In  this  section,  which  follows  [6),  the  step  transform  is  described 
mathematically.  Figure  45  shows  the  frequency  versus  time  diagrm  for  the 
received  linear  FM  signal. 
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g( t-uA) 


.  T 


uA ) ' 


(87) 


where  t  «  uA  is  the  time  of  arrival  of  the  received  signal,  B  is  its 
bandwidth,  and  T  is  its  duration.  Also  shown  in  the  diagram  is  a  segment 
of  the  sawtooth  reference  signal  used  to  Jeramp  the  input  signal.  The 
reference  segment  located  at  t  -  nA  is  given  by 


n*(  t-r.A) 


-5— ( t-nA) 
e  T 


(88) 


It  ha3  a  duration  of  T',  i.e., 

ml  mi 

-  —  +  nA  <  t  <  —  +  nA  . 
2  2 


After  performing  the  deramping  within  this  duration,  we  have 

^i(t-uA)2  t-nA) 2 

g(t-uA)gMt-nA)  ■  e  T  e  T  •  (89) 


Since  the  processing  is  performed  in  discrete  form, 

t  *  (nA  -  — )  +  kA,  0<k<  —  -1, 

2  A 


where  the  term  in  parenthesis  is  the  beginning  of  the  reference  segment,  k 
is  the  sampling  index,  and  A  is  the  sampling  interval.  The  discrete  form 
of  (89)  is 


g[(n+k-u)A  -  I-  l.g*(kA  - 

2  2 


(n+k-u) A 


Ill2  -iIB-(kA-  H)2 

2  e  T  2 


jH(u2A2+uAT* ) 

T 


j-T’I<n2A2-nAT* ) 
T 


,2r tBA  , 
j - (n-u)k 


2^8.2. 


un 


(90) 


This  output  is  then  Fourier  transformed,  with  respect  to  k,  to  extract  an 
image  with  a  resolution  corresponding  to  the  bandwidth  of  the  reference 
segment . 

The  interpretation  of  each  exponential  factor  in  (90!  is  as 
follows: 

i)  The  first  exponential  factor  is  independent  of  the  location  of  the 
sub-apeitr.re  reference  n,  and  is  solely  a  function  of  the  location 
of  the  input  signal  u,  therefore  this  term  will  not  be  affected  by 
the  first  or  subsequent  Fourier  transforms. 
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ii)  The  second  exponential  factor  is  an  undesirable  tent  3ince  it  does 
not  relate  -  to  the  input  signal  location,  u.  This  term  is  deter¬ 
ministic  anu  nust  bo  cancelled  out  before  the  second  Fourier  trans¬ 
form  is  taken. 

iii)  The  third  exponential  factor  is  a  linear  phase  term,  which  deter¬ 
mines  where  the  signal  spectrum  is  located  in  the  frequency  dimen¬ 
sion  after  the  first  DPT  is  taken  with  respect  to  k.  The  first  DFT 
can  be  regarded  as  a  filter  bank  which  segregates  the  spectral  com¬ 
ponents  in  the  sub-aperture.  The  spectral  locations  are  determined 
by  the  corresponding  locations  relative  to  the  sub-aperture  refer¬ 
ence. 


iv)  The  fourth  exponential  factor  is  a  phasor  whose  frequency  of 
rotation  depends  on  the  location  u.  In  order  to  recover  the 
precise  target  location,  another  Fourier  transform  can  be  applied. 


The  DPT  of  (90)  with  respect  to  k,  (ignoring  the  first  exponential 
factor)  is, 

1 2 


j^(n2A2-nAT  * ) 
X(r,n)  -  e  T 
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(91  ) 


whore  f{}  is  a  Fourier  transform  operator,  N  -  —  is  the  number  of 

A 

points  in  the  DFT,  r  is  the  frequency  variable  corresponding  to  the 
temporal  variable  k,  and  sinc( . )  is  a  cyclic  sinc-function  defined  as 
sinc(x)  -  sin(  ux)/sin(  ttx/n)  . 


The  location  of  the  sinc-function  depends  on  the  spatial  separation 
between  the  target  and  the  reference  ramp,  (n-u)A,  i.e.  the  target  is 
located  at  r  ■  r' 


r* 


NBA2 

T 


(n-u ) , 


(92) 


relative  to  the  sub-aoerture  reference  located  at  t  *  nA. 
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Assuming  tht  the  input  signal  is  sampled  at  the  Nyquist  rate,  we 
have  BA»1 .  If  the  signal  is  oversampled  (i.e.  BA<1),  as  is  normally  done 
in  practice,  some  of  the  Pourier  coefficients  in  r  will  be  small, 
containing  only  noise  and  leakages,  but  this  will  not  affect  our 
analysis.  Hence,  (92)  becomes 

.  NA,  > 
r  •  — (n-u) , 


T ' 

— (n-u ) . 


(93) 


Furthermore,  if  we  set  the  number  of  sub-apertures  that  can  be  fit  into  a 
full  aperture  to  arso  be  N,  i.e.  T/T'  ■  N,  then  (93)  can  be  simplified  to 


r ' 


(94) 


In  summary,  the  aperture  data  J  s  divided  into  N  sv.o-apertures  of  N  samples 
each,  which  gives  a  total  of  N2  points  for  a  full  aperture  data  set. 
Equating  (92)  and  (94)  gives  the  value  of  N 


1 

(FM  rate) A2 


(95) 


Since  adjacent  sub-aperture  references  are  separated  by  n^  cells, 
the  sub-apertures  can  be  indexed  by  a  new  parameter  £., defined  as  l  ”  n/N. 
Equation  (94)  then  becomes. 


r 

N 

Substituting  (95)  and  (96)  into  (91),  we  have 

jr  ■: 


'  J 


X  ( r ,  £ )  -  0 


(96) 


(97) 


Thus,  (95)  enaut's  that  ‘■ho  tnal  history  associated  with  a  point  target 
lies  on  a  lint  with  v  ’  c  c^pe  (in  the  r  versus  l  plane)  and  passing 
throug  .  the  1  axis  at  u/n,  as  denoted  in  (96).  This  arrangement  ma/tes  the 
processing  more  efficient  since  further  processing  can  proceed  along 
successive  diagonal  lines  with  unit  slope,  and  no  interpolation  is 
needed.  If  FFTs  axe  to  be  used  to  perform  DFTs,  N  must  '■«  a  power-of-two 
integer.  Unfortunately,  this  imposes  a  stringent  restriction  on  the 
signal  characteristic,  namely  the  FM  rate  and  the  sampling  interval,  upon 
which  a  step  transform  can  be  accurately  applied.  This  restriction  can  be 
partially  relieved  by  overlapping  the  sub-aperture  reference  ramps  as  will 
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be  described  in  Section  8.2.2. 


To  obtain  the  full  resolution,  the  sub-aperture  data  represented  by 
(91)  must  be  added  coherently  over  the  entire  aperture.  The  processing  i3 
two-dimensional  (across  aperture,  and  within  aperture,  'r').  The 

coherent  summation  is  performed  via  an  IDFT  with  respect  to  i(=n/N)  along 
the  diagonal  line  r”£,  i.e. 

N/2-1  -jtr(i2-l) 

Y(q)  -  l  [XU,D  e  ]  e  N  (98) 

i— N/2 


where  q  is  the  time  domain  equivalent  of  the  sub-aperture  number  'l'.  The 
exponential  factor  in  parentesis  cancels  the  first  exponential  term  in 
(97).  Before  the  second  DFT  is  taken  in  (99a),  X(r,£)  should  be  rotated 
over  r  by  half  a  cycle  (i.a.  N/2)  in  order  to  ensure  that  the  negative 
axis  in  r  precedes  the  positive  axis  in  a  linear  continuous  manner.  In 
practice,  this  can  be  achieved  by  embedding  an  appropriate  linear  phase 
term  in  the  reference  ramp.  If  we  substitute  from  (97)  into  (98), 

we  get 


Y(g)  - 


_  j2rru&  j2ffq£ 

N/2-1  N  N 

£  e  sinc[n]  e  , 

1— N/2  N 


sine  [a]  J  e  N 


_  j27rul  _j2irql 


e  N  , 


(99a) 


-  sine  [•£]  sinc[q-u]. 


sine  [— ]  sinc[q-uj. 
N 


(99b) 


Hence  after  both  DFTs  ( forward  and  inverse  DFTs  are  collectively  referred 
to  as  DFTs  since  they  do  not  bear  the  usual  connotations  of  the  Fourier 
transforms),  a  target  located  at  t*uA  is  recovered  as  a  compressed  pulse 
denoted  by  the  second  sine  function  in  (99b)  at  q-u.  The  first 
sinc-function  in  (99b)  is  an  envelope  modulation  function  inherited  from 
the  rectangular  window  In  the  first  DFT. 


As  mentioned  previously,  the  input  data  are  assumed  to  be 
adequately  sampled,  therefore  there  is  no  aliasing  effect  in  the  first 
DFT.  However,  as  described  in  (99a),  the  Nyquist  criterion  is  violated  in 
the  second  DFT.  The  width  of  the  mainlobe  of  the  envelope  function  is  2N, 
whereas  the  sampling  rate  is  only  N,  i.e.  the  DFT  length  is  only  N, 

Hence,  areas  of  the  mainlobe  corresponding  to  (uj>—  are  folded  back  into 
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unambiguous  region  of  the  DFT  resulting  in  heavy  degradation  of  the  radar 
image.  The  aliasing  effect  is  shown  in  Figure  46.  Three  point  targets  A, 
8  and  C  are  shown.  The  top  diagram  shows  that  A,  B  and  C  are  unambiguous¬ 
ly  resolved  as  if  there  is  no  aliasing  problem  in  the  second  DFT.  The 
bottom  diagram  shows  the  realistic  case  where-  the  function  is  circularly 
periodic,  with  period  q. 


Fig.  46  -  Compressed  image  after  coherent  summation  of  sob- 

aperture  images  along  a  diagonal  line  in  the  processing 
matrix.  Three  point  targets  A,  B  and  C  sure  shown, 

(a)  assuming  no  aliasing,  (b)  with  aliasing. 

8.2.2  Effects  of  Windowing  and  Sob-aperture  Overlap  on  the  Step  Transform 

In  order  to  avoid  the  aliasing  effect,  the  signal,  exp{-?^u^.}  in 

(99a),  must  be  adequately  sampled  in  'l',  and  the  point  target  response  of 
the  first  DFT,  i.e.  the  sinc-function  in  (99a)  must  be  replaced  by  one 
with  a  negligible  sidelobe  level. 

To  suppress  the  sidelobe  level  of  the  target  response  of  the  first 
DFT,  windowing  can  be  imposed  on  the  sub-aperture  data  before  the  DFT  is 
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taken.  In  practice,  a  Hamming  window  is  used  to  suppress  the  sidelobes  to 
less  than  -43dB.  In  such  a  case,  (91)  becomes 

j^.(n2A2-nAT') 

X(r,n)  •  e  T  e  T 

.2ttBA2, 

3 - ( n-u)k 

f{  [0.5.i-0.46cos(2iJs±Sji.)]  e  T  }, 

;££(n2A2-nAT*)  -^I^A2un  2 

»  e  e  w[r-— (n-u)],(100) 

T 

Nc  is  now  replaced  by  to  signify  the  effects  of  windowing  and  over¬ 
lapped  sub- apertures,  and  W(r)  is  the  Fourier  transform  of  the  Hamming 
window,  which  is  given  by 

W(r)  «  0.54D(r)  +  0.23[D(r-1 )+D(r+1 ) ], 


D(r)  -  e  N 


sin(nr) 


To  increase  the  sampling  rate  in  ’A',  the  spacing  between  adjacent 
sub-aperture  references  needs  to  be  smaller  than  N.  This  implies  that 
they  have  to  overlap.  Let  5  be  the  ratio  of  the  number  of  cells  overlap 
to  the  number  of  cells  in  a  sub-aperture.  The  sub-aperture  index  'l'  is 
redefined  as 


l  m - - - 

N^n-o 


(101  ) 


where  the  factor  N(1-£)  represents  the  sub-aperture  reference  spacing  in 
number  of  ceils. 

If  the  processing  efficiency  for  the  second  DPT  is  to  be  main¬ 
tained,  the  locus  of  the  envelope  function  of  X(r,&),  W(r,A),  given  by 
(100)  needs  to  lie  on  a  diagonal  line  parallel  to  r*A,  as  in  the  previous 
case,  otherwise  interpolation  would  be  required.  This  discrete  space 
arrangement  can  be  achieved  by  modifying  (95)  to 

2  T 

H  - - 

(1-OBA 


(102) 


(1-£)(fM  rate)A^ 


Substituting  (101)  and  (102)  into  (100),  we  have 
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jiri[<  1-0  2-1]  -jlHi 

X(r, 2)  -  e  e  N?  w[r-2+ - U - ], 

NCo-;r 


(103) 


Note  that  the  diagonality  property  (r-2  for  the  signal  position)  is 
maintained  in  the  argument  of  W[ . ] . 

We  can  perform  coherent  summation  along  the  diagonal  line  r»2  via 
an  IDFT  with  Hamming  window: 


N/2-1  2+0  5 

y  [o.54-0 . 46  cos  (2  n-————)  ]  X(2,2) 
2=»N/2  N;: 


-jir4[(1-5)l-1] 


e  ,  (104a) 


^2wui  ^2trq2 

-  w[ — - - ]  l  [0.54-0 . 46 cos  (2irAi2.iL)  J  e  N£  e  \  , 

»5d-5)  t 


w[— 

Nr  ( 1  — c  j 


W[q-u] , 


w[ - — - ]  W[q-u|, 

NC(1-5)' 


( 1 04b ) 


This  result  is  analogous  to  the  previous  case  without  windowing  and  sub¬ 
aperture  overlap,  as  described  by  (99b).  The  sinc-functions  are  now 
replaced  by  the  Fourier  tranuform  of  Hamming  window.  The  compressed  pulse 
is  represented  by  the  second  W-function  located  at  q*u.  The  first 
w-funcrtion  represents  an  envelope  due  to  the  Harming  window  in  the  first 
DPT .  The  W-function  offers  a  uuch  lower  sidelobe  level  (<-43dB)  than  the 
sinc-funct  on  (-13dE).  Hence  in  most  caes  the  w-function  can  be  regrarded 
as  sidelobe  free. 


A  minor  penalty  introduced  by  the  W-function  is  that  its  mainlcbe 
width  (null-to-null)  is  four  cells  wide,  whereas  that  of  the  sinc-f unctior. 
is  two  cells  wide.  The  mainlobe  would  cause  more  mainlobe  aliasing,  if  it 
were  not  for  the  overlapping  of  the  sub-apertures.  Assuming  the  sidelobe 
level  is  negligibly  low,  the  W-function  in  (104a)  limits  u  to  the  region 

-2  cells  <  — - -  <  2  cells,  recalling  that  the  Fourier  transform  of  the 

n5(i-C) 

Hamming  windowing  has  a  null-to-null  mainlobe  width  of  four  cells.  The  N- 

point  Fourier  transform  allows  <  u  <  ~£.  Equating  thes  inequalities 

2  2 

gives  a  sub-aperture  overlap  ratio,  5  "  0.75.  For  point  targets  lying 

outside  the  region,  i.e.  I  .  r'f  h ,  I  >  2,  the  W-function  attenuation 

( Nr ( I -c,  1  I 
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becomes  so  large  (>43dB)  that  aliasing  becomes  unimportant.  The  above 
aliasing  mechanism  is  illustrated  in  Figure  47. 


T  T 

(!■» 


Fig.  47  -  Aliasing  rachr-M  :-m  Atm  to  ttam  envelop  function 

of  the  Hassling  ■  jadow,  (a)  with  aliasing  for  £  <  0.75, 
(b)  aLnlana  condition  with  no  aliasing,  *  m  0.75. 


The  processing  efficiency  can  be  further  improved  by  minimizing  the 
sub-aperture  overlap,  and  hence  reducing  the  processing  of  duplicate 
data.  The  amount  of  overlap  can  be  reduced  by  allowing  aliasing  in  the 
portion  of  the  data  tht  will  be  thrown  away.  The  alias  free  region  is 
defined  as  the  width  of  the  separation  between  adjacent  sub-aperture 


references,  i.e., 
aliased  regions. 


-  iLiJ — <  u  <  .  Targets  falling  in  the 

namely,  -2Nr  ( 1  -  £ )  <  u  <  — LL  and  -£■— — —  <  u  < 

v  ~  2  2  ~ 


d-U, 


are  recovered  from  the  unaliased  regions  of  adjacent  diagonal  DFTs. 
Figure  48  illustrates  how  allowing  aliasing  in  the  throwaway  region  can 
reduce  the  amount  of  sub-aperture  overlay.  In  the  diagram,  the  envelope 
modulation  function  (the  first  terra  in  (104b))  is  partitioned  into  five 
regions  (A  to  E )  „  The  unambiguous  width  in  q  is  kept  at  .  Region  c 
represents  an  unaliased  region,  from  which  valid  processed  data  are 
derived.  The  other  regions  produce  duplicate  (invalid)  processed  data. 
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valid 


3n<i-C>-£ 


;  ■ 0.8 


Fig.  48  —  processing  efficiency  by  allowing  al jawing 

into  discarded  regions. 


It  is  desirable  to  let  region  A  to  fold  into  region  D,  and  region  E  to 
fold  into  region  B,  in  order  to  pack  all  regions  compactly  within  N 
cells.  Regions  B  and  D  now  contain  aliased  data  and  will  be  discarded. 
In  this  case,  as  shown  in  the  diagram,  the  optimum  value  of  is  0.6. 
Data  from  region  C  is  retained  as  valid  data.  The  valid  data  are 
amplitude  corrected  by  multiplying  by  the  inverse  of  the  envelope  function 

(i.e.  - - -  ).  The  process  is  repeated  for  all  diagonal  data,  and 

w( — 3—  ] 

nE(t-5l 

then  all  valid  data  are  concatenated  to  form  an  image. 

In  cases  where  0.2 <£<0.6  for  a  selected  value  of  N^,  the  data, 
which  lie  on  a  diagonal  after  the  first  set  of  FFTs,  are  undersampled  as 
just  described.  This  can  be  remedied  by  using  another  set  of  sub-aperture 
references  to  interpolate  the  original  diagonal  data  array.  The  method  is 
illustrated  in  Figure  49,  The  top  diagram  shows  the  signal  of  a  point 
target  with  the  original  sub-aperture  reference.  The  middle  diagram  shows 
the  signal  with  another  sub-aperture  reference.  The  new  sub-aperture 
reference  is  sc  designed  that  it  is  delayed  by  half  a  period  from  the 
original  sub-aperture  sawtooth,  and  has  a  negative  frequency  offset 
corresponding  to  half  a  resolution  cell.  This  offset  reference  provides 
an  additional  sampling  process  to  obtain  the  in-between  value  of  the 
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Pig.  49  -  Two  seta  of  reference  rteps  for  cases  with 

0.2  <  K  <  0.6. 


original  diagonal  data  array.  Phase  correction  is  then  imposed  on  both 
diagonal  arrays  as  described  in  (47).  Because  of  th»  additional  frequency 
offset  in  the  second  diagonal  data  array,  an  additional  phase  correction 
r  ttB  T  1  i 

factor,  expl-j - },  is  required.  This  is  done  to  maintain  phase 

T  2 

coherency  between  the  two  diagonal  arrays.  The  data  from  the  arrays  are 
interleaved  (as  shown  in  the  bottom  diagram)  before  the  second  FST  of 
length  2N  is  taken. 

In  cases  where  0  <  C  <  0.2,  one  has  to  choose  a  larger  value  of  N 
or  tolerate  a  fair  amount  of  degradation. 


3 


8.2.3  Conpater  Siwnlatlon  Result  a 

Figure  50  shows  the  result  of  a  computer  simulation  for  a  point 
target,  which  had  a  time -bandwidth  product  of  384,  and  which  was  processed 


SAMPLE  NUMBER 

Fig.  50  -  Ooagareaaed  pulse  of  a  point  target  processed  with 
step  transform. 


with  a  sub-operture  of  length  of  32  cel’s,  overlapped  with  the  other 
sub-apertures  by  20  cells  (62.5%).  The  prominent  main  peak  has  been 
normalized  to  OdB.  Artifacts  are  noticeable  but  are  mostly  at  least  -40dB 
down  from  the  main  peak.  These  artifacts  originate  from  the  sidelobe 
structure  of  the  first  w-f unction  given  in  (104b). 

The  humps  on  each  side  of  the  mainlobe  are  due  to  spectral  leakage 
from  the  incompletely  filled  sub-apertures  at  both  extreme  ends  of  the 
diagonal,  i.e.  the  data  starts  and  ends  part  way  through  those  sub¬ 
aperture.  These  partially  filled  sub-apertures  have  lower  resolution  than 
the  rest,  and  the  data  corresponding  to  the  target  is  spread  widely  over 
the  'r '-dimension  after  the  first  DFT  is  taken.  This  data  does  not  get 
properly  incorporated  in  the  Fourier  transform  along  the  diagonal,  but 
instead  it  contributes  to  the  integrated  sidelobe  ratio.  The  shapes  of 
the  humps  are  determined  by  the  window  used  in  the  diagonal  DFT.  The 


peaks  of  the  humps  are  located  at  approximately 


n£(1-5) 


from  either 


side  of  the  main  peak,  where  C  is  the  portion  of  sub-aperture  overlap. 
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8.2.4  Effects  c*f  the  Coble  Term  on  the  Step  Transform 

The  foregoing  description  of  the  step  transform  assumes  that  the 
input  signal  is  strictly  linear  FM.  It  is  thelinearity  (in  frequency 
versus  time)  of  a  linear  FM  signal  that  allows  the  step  transform,  which 
follows  the  trend  of  the  linear  FM  slope  in  a  step-wise  fashion,  to  be 
used  to  perform  pulse  compression.  The  introduction  of  cubic  and  higher 
order  phase  terms  destroys  the  desirable  linear  relationship,  and  results 
in  unequal  deramped  frequency  step  intervals,  and  other  abnormalities, 
such  as  defocussing  in  the  output  of  the  first  DFT.  In  this  section,  as 
was  done  for  the  matched  filter  approach,  we  will  consider  the  effect  of 
the  cubic  phase  term,  and  seek  ways  to  remedy  the  problem  it  causes. 

In  the  following  analysis,  we  assume  the  form  of  the  received  signal 

to  be 


}^.(t-uA)2  +  j2ire(t-uA)3 

g(t-uA)  *  e  T  ,  -  — +uA  <  t  <  — +uA, 

2  ~  “  2 

(105) 

where  t-uA  is  the  time  of  arrival  of  the  mid-point  of  the  received  signal, 
A  is  the  sampling  interval,  B  is  the  second- order  signal  bandwidth,  and  T 
is  the  second-order  signal  duration.  The  first  phase  factor  represents 
the  second-order  azimuth  signal  (i.e.,  a  linear  PM  signal).  The  second 
phase  factor  represents  the  third-order  azimuth  signal,  seen  at  extreme 
adar  squint  angles  [2].  The  coefficient  e  is  in  cycles/sec3. 

A  segment  of  the  sawtooth  reference,  centred  at  t-nA,  is  used  to 
deramp  the  input  signal.  The  reference  is  shown  in  Figure  (8.21),  and  is 
given  by. 


-J^(t-nA)2  (106) 

g*  ( t-nA )  -  e  T 

It  has  a  duration  of  T',  i.e., 

-  -'+  nA  <  t  <  nA.  (107) 

2  ~  ~  2 

After  multiplying  the  input  signal  by  the  reference  function,  we 

have. 


g(t-uA)  g*(t-nA) 

5^-(t-uA)  “+j2ire  (t-uA )  3  -J^(t-uA)2 

-  e  T  e  T  .  ( 108? 

Because  the  processing  is  performed  in  discrete  form,  t  must  be 
replaced  by 

m  •  ml 

t  -  (nA  -  — )  +  kA,  0  <_  k  £  N^-1  ,  Nr  -  — , 

2  A 


(109) 
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where  the  terra  in  parenthesis  is  the  beginning  of  the  reference  segment, 
and  k  is  the  sampling  index  within  the  sub-aperture.  The  discrete  form  of 
(108)  is 


g[(n+k-u)A  -  — ]  g*(kA  -  —  i 
2  2 

;£i(u2A2+uAT* )  j^.(n2A2-nAT' )  n-u ) k 

-  e  T  e  T  e  T 

2 

un  j2ire[(n+k-u)A  -  -  l3.  (110) 

h  T  e  2 

The  first  factor  in  (110)  is  solely  dependent  on  u  and  is  of  no 
importance  to  the  image  processing.  It  will  therefore  be  dropped.  with 
this  change  (110)  becomes 

g[(n+k-u)A- — j  g*(kA-— ) 

,xB,  2,2  .2tr3A2,  ,,  ,2wBA2  r,,  ,  .  T '  .  .13 

j—  (n  A  -nAT')  j - (n-u)k  - j-"..— —  un  j2wel (kA-uA- — )+nAJ 

■  e  T  eT  e  T  e  2 

(111) 


The  last  exponential  factor  in  (111)  is  solely  due  to  the  cubic 
phase  term.  The  k-dependency  in  this  factor  signifies  that  the  error 
caused  by  the  cubic  term  is  proportional  to  the  duration  of  a 
sub- aperture.  In  practice,  the  sub-aperture  length  is  much  shorter  than 
the  full  aperture,  and  the  variation  in  the  phase  factor  as  k  changes  from 
one  end  of  the  practice,  the  sub-aperture  length  is  much  shorter  than  the 
full  aperture,  and  the  variation  in  the  phase  factor  as  k  changes  from  one 
end  of  the  sub-aperture  to  the  other  is  vary  small.  An  expansion  of  thi-* 
factor  gives 


j2ue[(kA-uA-— )+nA]3 
e  2 


j2irt(kA-uA' 

e 


3 


j2it[  3c  (kA-uA-—  )2nA] 
e  2 


j2ire(nA)3  j2n  [  3  e  (kA-uA-^— )  ( nA )  2  j 
e  e  2  .(112) 


The  significance  of  these  phase  terms  is  depicted  in  Figure  51,  and  ij 
explained  as  follows; 

i)  The  first  exponential  factor  is  independent  of  n.  It  represents 
the  phase  deviation  due  to  the  cubic  term  of  the  input  signal 
within  the  sub-aperture  located  nearest  to  the  centre  of  the 
full  aperture.  This  term  is  extremely  small  and  is  of  the  order 
of  less  than  0.5*  for  hignly-squinted  (e.g.  6*  from  flight 

direction),  high  resolution  {e.g.  1  metre),  and  long  slant  range 
(e.g.  100  km)  cases  operating  at  x-band. 
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ii)  The  second  exponential  factor,  which  is  a  function  of  n, 
represents  the  additional  phase  error  within  a  deraraped 
sub-aperture.  This  phase  error  is  rather  small  under  most 
stringent  circumstances,  it  is  less  than  10*  for  sub-apertures 
located  at  the  ends  of  the  full  aperture. 

iii)  The  third  exponential  factor  depends  solely  on  n.  ltdescri.bes 
the  global  phase  error  along  the  input  data.  Its  magnitude  may 
be  rather  significant  and  it  must  be  cancelled  before  the  second 
DPT  is  performed. 

iv)  The  fourth  exponential  factor,  depandent  on  k,  u,  and  n,  is  the 
amount  of  frequency  c-ffset  from  the  normal  decamped  CW  signal. 
This  frequency  offset  shows  up  as  a  shift,  in  the  spectrum 
produced  by  the  first  DPT. 

In  summary,  the  first  tvt  exponential  factors  describe  the  local 
residual  phase  errors  within  a  sub-aperture  Their  -Magnitudes  are 
normally  too  small  to  affect  the  final  result.  The  lajh  two  exponential 
factors  describe  the  global  residual  pht3e  error  and  the  frequency  offset 
respectively.  Tueir  effect  on  the  final  result  is  worth  further 
consideration. 


> 


FM  SIGNAL 

Pig.  SI  -  Deramping  a  non-linear  PM  signal  with  coble  phase.  All 
errors  are  shown  for  the  case  in  which  an  ideal  sawtooth 
reference  is  applied. 


With  only  the  last  two  phase  factors  of  the  expansion  included,  the 
decamped  signal  of  (6.30)  can  be  written  as 
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r  T' i  T ' 

g  j  ( n+k -u)  A> - I  a*(kA- - ) 

2  2 


}^(n2A2-nAT-)  }ilR. 
e  T  e  T 


*RA2,  ,2itBA2  .  „,3 

- -(n-u)k  -  j— — -  ■  un  }2TTe(nA) 

T  e  T  e 


j2it[3e(kA-uA-H)  (nu)2] 


M13) 


The  image  is  produced  from  <  1  » 3 )  by  using  the  modified  step 
transform.  First  a  Hamming  window  is  imposed  on  the  sub-aperture  in  order 
to  suppress  sidelobe  leakage,  and  a  DFT  is  taken  with  respect  to  k?  the 
result  13 

(n2A2-nAT ' )  j2ne(nA)3  -j6ne(uA-t— )  (nA) 2  -  j-*'11—  un 

X(rr,n)  »■  e  T  e  e  ^  e  ^ 


N-1  k+0.5 

l  [0. 54-0. 46  003 ( 2n - )  ] 

k-0 


[e  T 


-(n-u)k  j2s3ekA(nA) ‘ 


]  e  , 


^5.(n2A2-nAT' )  j2*e(nA)3  -j6*E(uA-»2— )  (nA)2  un 

k  ^  o  A  2  A  T 


w[r-( 


NrBA^(n-u) 


+3N  EA3n2)]. 


(114) 


where  w(r)  is  the  Fourier  trarsform  of  the  Hamming  window  and  is  given  by 
W(r)  -  0. 54D ( r )-0 .23 [0( r-1 )+0( r+1 > ) ,  (115) 


D(r )  -  e 


■Jjr  3in[nr] 

,  irr , 
sinl — ] 


Substituting  (101),  (102)  and  T'  «-  N^A  into  (113)  yields 


-j6TiCN|(1-5)2A3uZ2  jnRU) 


X(  r ,  f.)  -  e 


Nr  r  3  7  3^ 

w[r-Jl-3£Hf ( 1  -U“A  i'-1 


(116) 


NU1-S) 
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where  6(1)  -  2e [Nr { 1 -i ) &] 3 23  +  [ 1-e- sen3 (1 ) 2A3 ] l2  l.  We  can  see  that 

5  n  n  ^ 

the  location  of  the  spectrum  is  off  by  [3cN^(1-5)‘,4  r]  cells  from  the 
normal  diagonal  position.  If  uncompensated  thi„  iOaitional  shift  will 

result  in  amplitude  tapering  of  the  sigr  1  on  both  ends  cf  the  diagonal; 
moreover,  the  energy  associated  with  u.ie  diagonal  can  curve  into  the 
adjacent  diagonals.  This  intermixing  of  diagonals  will  introduce  image 
blurring  into  the  DPT  of  the  adjacent  diagonals. 

If  this  amount  of  degradation  is  intolerable,  then  the  curvature 
can  be  corrected  by  interpolation  along  the  r~d^mension  after  the  first 
DPT. 


The  W-function  in  (116)  represents  a  narrow  bandpass  filter  moving 

3  2  4  2 

along  a  locus  described  by  r  *  l  *  3eN^(!-C)  A  £  ,  in  the  r  versus  £ 

plane.  To  regain  the  original  resolution,  a  coherent  summation  (DFT)  over 
r  needs  to  be  performea  along  this  locus. 


The  first  exponential  factor  in  (116) is  dependent  on  the  location 
of  the  target,  u,  and  the  sub- aperture  location,  2.  The  maximum  phase 
value  of  this  term  occurs  at  the  above  extreme  values 


i.  i 
i  u  I 


NC(1-U  ,  NC 

—  CUiU  U|  -  - - 

2  '  I  2 


9max  ■  0.75tteN5a|(1-?!3, 

0.75xeT3 
“  - - 


(117) 


The  maximum  tolerable  ^niax  value  is  0.75tt  for  heavily  weighted  DFT 
(3ee  (120)  bel'.iw).  Therefore 

>  £T3 .  (1 18 ) 

Substituting  the  criterion  in  (118)  into  the  cubic  phase  at  the  end  of  a 
full  aperture  as  given  in  (105),  gives 

9cubic  “  )  ' 

dm 


TtNr 

<  — (119) 

4 

In  other  words,  if  the  cubic  phase  at  the  end  of  the  full  aperture  is 
xN^/4  or  less,  then  the  first  exponential  factor  in  (116)  is 
insignificant,  and  will  be  dropped  henceforth.  For  example  if  Nr  »  32, 
the  maximum  tolerable  cubic  phase  at  the  end  of  the  full  aperture  is  1440; 
a  phase  error  which  would  otherwise  be  intolerable  if  no  sub-aperture 
compensation  is  used. 
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Before  proceeding  with  the  coherent  sub-aperture  summation,  namely, 
a  second  DPT,  the  second  and  third  exponential  factors  in  (116)  should  be 
cancelled.  When  this  i3  done  the  second  DPT  with  Hamming  window  becomes. 


N?/2~1 


j2rrq£ 

J  *1 


(q)  -  7.  X(r,f)  [0.54-0.46cos(2itH^ll)]  e  , 

Nr/2 


wf - - —  ]  r  [o  .  54-0 . 46cos  f  2 

•NCd-e>  i  ' 

along 

r-t+3eNf  ( 1  - ^  )  2A3 ii2 


.2nul  ,2itqS, 
■1 -  3 — — 

Nr  Nr 


„[ - ^ - ]  w[q-u  ] . 

n4(i-£) 


(120) 


After  both  DFTs ,  a  target  located  at  t  *"  uA  is  recovered  as  a  com¬ 
pressed  pulse,  w[ q-u ] .  The  first  W-function  in  (120)  represents  an  unue- 
sired  amplitude  modulation  function  inherited  from  the  Hamming  weighting 
in  the  first  DfT,  and  should  be  corrected. 


8.2.5  Range  Curvature  Compensation 

Range  curvature  can  be  compensated  in  a  similar  fashion  to  the 
matched  filtering  case  described  in  Section  7.1.  In  the  step  transform, 
the  range  curvature  compensation  should  take  place  after  the  first,  set  of 
DPTa.  Referring  back  to  Figure  42,  each  row  of  data  indexed  in  l  have  the 
same  frequency  characteristic  and  are  acquired  through  the  same  portion  of 
the  antenna  pattern.  Thus  each  row  of  data  in  the  prot  -'ssing  matrix 
Buffers  the  same  amount  of  curvature,  whereaa  tne  data  in  a  column,  which 
are  each  acquired  from  different  angles,  suffer  different  amounts  of 
curvature,  To  perform  the  interpolation  a  number,  based  on  the  length  of 
the  interpolator,  of  processing  matrices  from  consecutive  ranges  are 
collected  (see  Figure  52).  Interpolation  is  then  performed  across  the 
matrices  on  a  set  of  corresponding  dcua  (f,r),  one  from  each  processing 
matrix. 


‘J  .  COHK/TATIOWAb  RBQUrJtEHiiHTS 

In  this  section,  we  attempt  to  gi’e  an  assessment  of  the 
computational  requirements  for  the  foregoing  SAR  data  processing 
algorithms.  SAR  data  processing  ia  often  characterized  by  its  voluminous 
data  size,  and  the  stringent  requirement  for  antlimotic  operation?.  Inus 
the  design  of  a  SAR  processor  should  lx:  optimized  in  terms  of  arithmetic 
complexity  and  control -function  complexity.  The  predominant  factors  for 
consideration  are  the  throughput  rate,  the  poise-compress  ion  ratio,  ai,u 
the  number  of  looks.  A  deairabLe  algor i th;n  for  SAR  data  processing  should 
be  selected  by  making  the  necessary  tradeoffs  between  these  factors. 
Because  of  a  large  number  of  parameters,  *  he  cninputat  iona 1  rates  of 
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m*1  m  m-1 


CONSECUTIVE 

RANGE 

CELLS. 


rnoocoon'iu  mmnlCtS 

OF  CONSECUTIVE 
RANGE  CELLS 


J 


Fly.  52.  -  Bunge  auvatace  ca^pensatloa  in  step  transform. 


generalized  optimal  processing  schemes  are  not  easily  compared.  In 
general,  the  matched  filtering  approach  requires  less  control-function 
complexity  and  provides  more  flexibility  and  exactness  in  tailoring  the 
matched  filter  to  fit  the  data.  On  t'  «  other  hind,  the  deramping 
technique  demands  less  memory  and  can  proceed  in  a  continuous  rather  than 
batch  fastiior,  and  is  morj  cuited  to  real-time  processing.  In  the 

present  asnesament,  v.  re3s  the  computational  requirements  in  terms  of 
numbers  *  complex  multi*.* ^cations  and  additions  (including  subtractions). 
In  the  hope  that  the  quantity,  when  considered  with  other  factors,  will 
help  in  the  refinement  of  the  overall  processing  configuration.  It  is 
assumed  in  the  following  that  the  matched  filters,  or  the  deramp 
references,  and  the  phase-magnitude  correction  factors  are  pro-computed, 
and  that  the  effort  required  to  synthesize  them  is  small  compared  with  the 
actual  processing  of  the  data. 

Unless  otherwise  redefined,  the  symbols  used  in  the  following 
evaluations  are  the  same  as  those  used  in  their  respective  descriptions 
given  earlier.  In  order  to  use  the  FPT  algorithm,  the  variables  N.j , 
tls  and  N  need  to  power-of -two  integers. 

1  Cca.Aitf  i.  '  'Requirement. a  tor  Single  book  Case*  without  Range 
Curve*.  w.  ubic  Phase  Term  Ccanj reum s t  i 

"ne  computai  icnal  requirements  for  the  three  pui.se  compression 
techn*.  as  i  eva  uateri,  Tt  is  assuina d  ti  at  the  complications  due  to 
range  tui  vuuuo  and  th*1  cu'-K  pi-  Lena  ai.e  r.-«»g  i  igible. 
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Computational  requirement  for  the  matched  filtering 
approach  (via  fast  convolution) 


For  complex  multiplications 


-Forward  FFT  on  a  set  of  My  input  data 


-Matched  filtering  with  windowing  included 


-Inverse  FFT  on  My  matched  filtered  data 


llog2NT 


-Total  number  of  complex  multiplications 


N«y  log2“^T^^T 


For  Complex  Additions 


-Forward  PFT  on  My  input  data 


Mr 

—  --1  og  2 


-Inverse  PFT  on  t*y  matched  filter  data 


Log  2  My 


-Total  number  of  complex  additions 
-Matched  filter  length  in  time  domain 
-Valid  output  data  length 


N«ylog2My 


-Number  of  complex  multiplications  per 
valid  output  point 


log2NT+1 


-Number  of  complex  additions  per  valid 
output  point 


log2NT 


b Computatio  nal  requi  rement  for  the  step  transform 

Let  !_  be  the  sub-aperture  overlap  and  Nf  be  the  number  of 
uub  apertures,  so  that  N’f  ■ - - - - - , 
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For  cages  with  C>0.6 

For  complex  multiplications 

-Deramp  input  signal  with  Nj:  sub-apertures 
and  use  data  points  for  each  sub- 
aperture 

-First  FFT  for  one  sub-aperture 


-Do  the  above  for  Nj;  sets  of  3ub-aperture 
FPTs  which  results  in  a  NjrXN^  matrix 
containing  N^  complece  diagonals 

-Number  of  multiplications  to  form  one 
complete  diagonal 


-ilog2Nc+Nc 


—Phase  our  tec  i.  ion  eu>d  windowing  on 
one  diagonal 


-Second  FFT  (along  diagonal) 


Nr 

— — log2Njr 


-Total  number  of  complex  multiplications  Hrlog2F^>2N^ 

(not  including  the  final  amplitude 
correction) 

For  complex  additions 

-Number  of  additions  in  the  first  and 
second  FFTs  that  are  required  to  form  N^log2N^ 

the  output  from  one  diagonal  - 

-Number  of  valid  data  produced  from  (l-^)N^ 

one  FFT  along  the  diagonal 


-Number  of  complex  multiplications 
per  ■'ralid  output  point 


log2N£+2 


-Nuitiber  of  complex  additions  pe.; 
valid  output  point 


log2N£ 

~~-T 
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For  cases  with  0.2<£<0.6 


These  cases  require  two  sets  of  sub-aperture  references. 
Section  8.2.2  for  details. 

For  complex  multiplications 

o 

-Deramp  input  signal  with  2Njr  2N^ 

sub-apertures  and  use  data  points 
for  each  sub-aperture. 


-First  FFT  for  one  sub-aperture 


-Do  the  above  for  2N^  sets  of 
sub-aperture  PFTS  which  results  in 
2N^  X  matrix  containing 
complete  diagonals 

-Number  of  nultlplications  to  form 
one  complete  diagonal 

-Phase  correction  and  windowing  on 
one  diagonal 

-Second  FFT  (along  diagonal) 


N£log2N£+2N£ 


2N*r 

2Nr 

—  s.log22Njr 
■Njr  ( 1+log2N^) 


-Total  number  of  complex  multiplications 
(not  including  the  final  amplitude 
correction) 

For  complex  additions 


-Number  of  additions  in  the  first  and 
second  PFTs  that  are  required  to  form 
the  output  frcm  one  diagonal 


2N£l092+SN£ 


N£log22N^ 

■2Nj- log2N£+N£ 


-Number  of  valid  data  produced  (1-£)Nr 

from  one  FFT  along  the  diagonal 


-Number  of  complex  nuitipl ications 
per  valid  output  point 


— 

2N^log2Njc+5 


see 
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Number  of  complex  additions  per 
valid  output  point 


21og2N£+1 


c)  Computational  requirements  for  the  spectral  analysis  method 


In  the  present  treatment,  we  assume  that  the  guard  band,  0,  and  the 
DFT  length,  N,  depend  on  the  required  resolution,  and  that  they  are  suf¬ 
ficiently  large  that  the  number  of  processing  columns  within  a  parallelo¬ 
gram  is  not  less  than  two  (see  for  example.  Figure  38).  The  number  of 
processing  columns  required  for  each  parallellograra  can  be  calculated,  by 
using  (82),  as 


L  -  INTEGER  [ 


N 


a 


]. 


•  INTEGER  [ 


+  1 


1-6' 


oNk 


Let  Np  be  the  recorded  data  length  before  presumming,  therefore 
B2 

N  «  BT  -  -  . 

H  k 


For  complex  multiplications 

-Deramp  input  signal  (including  windowing 
to  form  L  processing  columns  within  a  LN 

parallelogram ) 

-FFTs  for  L  processing  columns  L^lc^N 


-Total  number  of  complex  multiplications 
for  one  parallelogram. 


For  complex  additions 

-Number  of  additions  in  the  L  FFTs 


i?-llog2N+NL 


LN. 

— log2N 
2 


-Number  of  valid  data  given  by 
one  parallelogram 
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-Number  of  complex  multiplications 
per  valid  data  output 


-Number  of  complex  additions  per 
valid  data  output 


9.1.1  Examples  for  Cases  without  Range  Curvature  and  Cubic 
Phase  Tern  Co^ensatioc 


The  following  examples  are  intended  to  provide  some  insight  into 
the  computational  requirements  of  the  processingalgorithms .  Two  cases  are 
given,  both  have  a  prasum  factor  of  8. 


Examples  for  typical  resolution  (3.5m) 

Parameters: 

-Raw  data  length 

-Data  length  after  presum, 

-Slant  range 

—Matched  filter  length,  M 
( >Time-bandwidth  product  in  cycles) 

a )  Matched  filtering  approach 

Number  of  complex  multiplications 
per  valid  data  point 


Number  of  complex  additions 


Ml 


16,384 
2,048 
100  km 
557 


(log22048)+1 


2048 


16.5 

log22048 

1 -Hi- 

2048 


15.1 


b)  Step  transform  (assuming  (8.21)  is  satisfied) 

The  computational  requirement  depends  heavily  on  the  FM  rate  and 
the  sampling  interval  as  given  in  (0.21),.  or  equivalently,  the 
oversampling  factor.  According  to  the  Nyquist  criterion,  BA<1 ,  we  have 

N,;  > 


TE> 

1-5 
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a nd 


<  1  - 


TB 

4 


where  TB  is  the  time-bandwidth  product.  Setting  TB  ■  557,  £  ■  0.6,  we 
have, 

N,>  HE,  37.3, 

*  /  0.4 

we  choose  *  64,  and 

5  <  3  -  -”I2 

~  ( 64  ) z 

<  0.864. 

which  gives  the  number  of  overlapped  points  ■  27.6  *  27,  £  is  recalculated 
as  0.844 . 

To  calculate  the  computational  requirement  with  £  •  0.844,  we  have 


-Number  of  complex  multiplications 
per  valid  data  point 


log; 


64+2 


1-0.844 


52 


-Number  of  complex  additions  per 
valid  data  point 


log; 6 4 
1-0.844 


-  _39 

c)  The  spectral  analysis  method 

Since  N  -  557,  we  choose  N  ”  1024  ,  this  makes  the  resultant 

resolution  proportionally  finer. 

The  number  of  processing  columns  is 

T.  -  INTEGER  [  - - - +  l],  where  8“0.4,  I;*''.. 4, 

1.-0. 4  .  0^4*1024 
16384 


-  INTEGER  [2.7], 
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-Number  of  complex  imiltiplications 
per  valid  data  point 


—logo  1 024+2 
2 


-Number  of  complex  additions  per 
valid  data  point 


-dog?1024 

2 


Examples  for  fine  resolution  (1.75  m) 


Parameters : 

-Raw  data  length 

-Data  length  after  presum 

-Slant  range 


16,384 
2,048 
100  km 


«.  i.  .U  —  J  £1  1i _ 1 - 1-  %- 

"nacutcu  li-itu.  tcn^Lii 


1,116 


Matched  filtering  approach 


-Number  of  complex  multiplications 
per  valid  data  point 


(log22048)+1 
.  1116 


-Number  of  complex  additions 
per  valid  data  point 


log22048 

1116 


b)  Step  transform  (assuming  (102)  is  satisfied) 


2  TB 
Nr  >  - 

4  “  1-5 


>  52.8 


choose  N(r  -  64 


C  <  1 


fl 
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£  0.728 

which  gives  the  number  of  overlapped  points  ■  2.33  *  23  .  C  is  recal¬ 
culated  as  0.719. 

To  calculate  the  computational  requirement  with  £  -  0.719,  we  have 


-Number  of  complex  multiplications 
per  valid  data  point 

-Number  of  complex  additions 
per  valid  data  point 


log264+2 

1-0.719 

29 

log264 

1-0.719 


22 


c)  Spectral  analysis  method 

N  -  1116 
choose  N  ■  2048 

The  number  of  processing  columns  is 

L  -  INTEGER  [ - 1 - +  ll,  with  6  -  0.4,  a  -  0.4. 

1-0.4-.?-?*2048 

16384 


-  INTEGER  [2.8  t 

-  2 

-Number  of  complex  nu1 .iplications  per 
valid  data  point 


— log22048+2 


-  13 

-Number  of  complex  additions  per  2. 

— xog22048 

valid  data  point  2 

-  _1_1_ 

In  practice,  the  matched  filtering  method  acid  the  step  transform 
reauire  the  raw  data  be  presummed  before  being  processed.  The  presumming 
can  be  accomplished  by  using  a  two-stage  low- pass  filter  with  15  coeffi¬ 
cients  each.  Therefore,  an  additional  30  multiplications  and  additions 
per  input  data  point  should  be  added  on  top  of  the  computation 
requirements.  In  the  case  of  the  spectral  analysis  method,  the  presumming 
process  is  done  by  selectively  derampinq  and  windowing  the  desired 
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processing  coluiim,  thus  no  additional  presumming  is  required.  However, 
the  output  product  ia  not  as  well  controlled,  wich  respect  to  aliased 
energy,  as  the  cases  in  which  a  presumming  filter  is  used. 

9.2  Computational  Requirements  for  Single  Look  Cases 
with  Range  Curvature  Compensation 

The  above  computational  evaluations  are  now  extended  to  include 
cases  with  range  curvature.  The  spectral  analysis  method  is  not  included 
since  it  cannot  accommodate  range  curvature. 

a )  Computational  requirements  for  the  matched  fi3tering 

approach  (via  fast  convolutional)  with  range  cu rvature 
compensation 

The  range  curvature  compensation  is  accomplished  by  interpolation, 
in  the  range  dimension,  along  the  curvature.  A  four-point  interpolator  is 
assumed  here.  The  procedure  has  been  described  in  Section  6.2 

For  complex  multiplications 

-4-paint  interpolator  4NT 

-Other  computations  NTlog2N*i>+NT 


-Total  number  of  complex  nultiplications 

For  complex  additions 
-4-point  interpolator 
-Other  computations 


N^log^N'ptSN'p 


Nijilog  2N>ji 


-Total  number  of  complex  additions 


NTlog2N+3NT 


-Valid  output  data  length  N^i-M 


-Number  of  complex  multiplications  per 
valid  data  output  point 


log2NT+5 


Nr 
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-Number  of  complex  additions  per 
valid  data  output  point 


log2  Nij+3 


b)  Computational  requirements  for  the  step  transform 
with  range  curvature  compensation 

The  range  curvature  compensation  is  accomplished  by  interpolation 
of  the  corresponding  elements  in  the  processing  matrices  from  different 
ranges,  as  described  in  Section  8.6.  A  four-point  interpolator  is  assumed 
here  also. 


For  cases  with  5  >  0.6 
For  complex  multiplications 

4N? 

Nt-log2N£+2N£ 


-Interpolation  to  form  one  complete  diagonal 
-Other  computations 


-Total  number  of  complex  multiplications 


N^log2N^+6N^ 


Por  complex  additions 


-Interpolation  to  form  one  aomplete  diagonal  3N^ 


-Other  computations 


Nj-log2N£ 


-Total  number  of  complex  additions 


N^log2Nf +3Nr 


-Number  of  valid  data  given  by  one  (1-£)Nr 

diagonal  FFT 


-Number  of  complex  nultiplications 
per  valid  data  output  point 


log2N^  +6 


-Number  of  complex  additions 
per  valid  data  output  point 
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For  cases  with  0.2  <  5  <  0.6 


For  complex  multiplications 

-Interpolation  to  form  one  complete 
diagonal 

-Other  computations 


8N,r 

2N£log2N£+5N£ 


-Total  number  of  complex  multiplications 


2N ^ log2Nr  + 1 3N  jt 


For  complex  additions 

-Interpolation  to  form  one  complete  diagonal  6N^ 

-Other  computations  2N£log2Nt;+N5 


-Total  number  of  complex  additions 


2Njrlog2Nr+7N£ 


-Number  of  valid  data  given  by  one  (1-5  JN^- 

diagonal  FFT 


-Number  of  complex  nultiplications 
per  valid  data  output  point 


21og2Nr+U 

W; 


-Number  of  complex  additions  per 
valid  data  output  point 


21og2Nj-  +7 

1-5 


9.3  Computational  Requirements  for  the  Coherent  Multi-look 

Cases  with  Range  Curvature  and  Cubic  Phase  Compensations 

The  computational  requirements  for  the  severe  combination  of  range 
curvature,  cubic  phase,  and  high  resolution  are  encompassed  in  the 
f ol lowing . 


I 


Computational  requirements  for  the  multi-look  matched 
filtering  approach  (via  fast  convolution) 


No  additional  computation  is  required  for  the  cubic  phase  compensa¬ 
tion,  since  the  cubic  phass  tent  can  be  pre-caiculated  and  embedded  i-.  the 
matched  filter  when  it  is  synthesized.  If  sub-aperture  summation  is  per¬ 
formed  in  the  time-domain,  the  inverse  Fourier  transform  of  each  sub- 
aperture  has  to  be  the  3ame  length  as  the  overall  original  data  length 

(Nt).  Thi3  is  accomplished  by  appending  zeroes  to  the  unoccupied  por¬ 

tion  of  a  sub-aperture  spectrum.  If  sub-aperture  summation  is  performed 
in  the  frequency  domain  during  compression,  the  sui-aperture  spectra  are 
concatenated  (except  in  the  transition  bands)  to  the  original  bandwidth, 
before  the  inverse  Fourier  transform  is  taken.  Thus  the  amount  oi 

computation  is  very  much  the  same  as  that  for  the  single  look  case.  In 

the  following  analysis,  post-compression  sub-aperture  summation  in  the 
time  domain  is  assumed. 


For  complex  multiplications 

-Forward  FFT  an  a  set  of  NT  input  data 

-Interpolation  for  all  sub-apertures 
(with  a  4-point  interpolator) 


-Matched  filtering  with  extraction  filter 
and  windowing  for  L  sub-apertures 

-L  inverse  FFTs  with  NT  date 


-Total  number  of  complex  multiplications 


?1 092^1 


i-^Xlog^N'T 


Nt(L+1 ) 


log2N'j' 


For  complex  additions 


-Forward  FFT  on  a  set  of  input  data 


-Interpolation  for  all  sub-apertures 


-i,  inverse  FFTS  with  N  data  each 


•Coherent  iumeution  of  I.  sub-apartuves 
in  the  '.oa'.'il  domain 


-'-^log  ?  N-v 
-\  *■  * 


-log  ^  N'r 


(L-1  )  N-, 


m 
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-Total  number  of  complex  additions  Nij!L+1  ) 

- log^N^. 

*{  L+2  )NT 


-Valid  output  data  length 


Nt-M 


-Number  of  complex  multiplications  per 


valid  output  point 


-Number  of  oomplex  additions  per 


^-log2NT+I.+2 


b)  Computational  requirements  for  the  step  transform  with 
range  curvature  and  cubic  phase 


If  the  cubic  phase  terra  is  significant  compared  to  the  quadratic 
phase  term,  then  the  data  array  for  the  second  DFT  will  no  longer  lie 
along  the  diagonal  of  the  array.  This  phenomenon  has  been  described  in 
Section  8.6.  To  remedy  this,  interpolation  is  used  to  straighten  out  the 
curvature  of  the  diagonal,  before  the  computation  of  the  second  DFT.  In 
the  following  analysis,  we  assume  that  all  the  diagonal  data  are 
interpolated  using  a  four-point  interpolator. 

For  cases  with  C  >  0.6 

For  complex  multiplications 

-Interpolation  for  the  diagonal  data  4N^ 

-Other  calculations  (including  range  N^log2N^+6N^ 

curvature  computation) 


-Total  number  of  oomplex  nultiplications  Nf log2N^ +1 0N^ 


For  complex  additions 


3N£ 


-Interpolation  for  the  diagonal  data 


-Other  calculations  (including  range 


N^log2N^+3N^ 


curvature  compensation) 

-Total  number  of  complex  additions 

-Number  of  valid  data  given  by  one 
diagonal  FPT 


Number  of  complex  multiplications 
per  valid  data  output  point 

Number  of  complex  additions 
per  valid  data  output  point 

For  canea  with  0.2<5<0 .6^ 

For  complex  multiplications 

-Interpolation  for  the  diagonal  data 

-Other  calculation  (including  range 
curvature  computation) 

-Total  number  of  complex  multiplications 
For  complex  additions 
-Interpolation  for  the  diagonal  data 

-Other  calculations  ( including  range 
curvature  compensation) 

-Total  number  of  complex  additions 

-Number  of  valid  data  given  by 
one  di  agonal  FFT 


-Number  of  complex  multiplications 
per  valid  data  output  point 


-Number  of  complex  additions  per 
valid  data  output  point 


N^log2^r+6N^ 


m-C)n5 


log2NF-*-10 

i-e 


log2N|r+6 

1-5 


BN? 

N^log22N|+12N5 

N^log22N|+20N^ 

6N^ 

N£log22N£+6N£ 

N^log22Nr+12N? 


<1-£)NC 


log2  2n| -t-20 
1-5 


log22NrH-12 

1-5 
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10.  SUMMARY 

A  general  theory  of  coherent  sub-aperture  processing  for  SAR,  has 
been  presented.  The  problems  encountered  in  applying  a  one-dimensional 
matched  filter  on  SAR  azimuth  data  have  been  reviewed  and  cnaracterized. 
A  piecewise  solution  for  the  one-dimensional  matched  filter  has  been  form¬ 
ulated  and  it  has  been  shown  how  this  solution  can  be  used  to  combat  the 
processing  problems  without  resorting  to  memory  intensive  two-dimensional 
filtering.  The  solution  has  been  thoroughly  tested  by  computer  simula¬ 
tions,  and  an  error  analysis  of  the  effects  o‘  the  phase  and  amplitude 
distortions,  caused  by  the  piecewise  approximation,  hats  been  given.  In 
addition,  the  basic  deramping  techniques  have  been  overviewed,  and  the 
step  transform  technique  discussed  in  detail.  In  particular,  the  step 
transform  method  has  been  extended  to  include  compensation  for  phase 
errors  and  range  curvature.  Finally,  the  computational  requirements  for 
different  cases  and  processing  techniques  have  been  evaluated. 
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APPENDIX  A 


Polynomial  Truncations  of  the  Taylor  Expansion  of  the 
Radar  Received  Signal 

To  illustrate  that  the  quartic  (and  higher  order)  phase  terns  are 
insignificant,  we  go  through  the  following  numerical  example. 


We  assume  a  high  squint  angle,  high  resolution  case,  which 
represents  the  worst  possible  operating  conditions: 


Resolution  (p  )  - 

wavelength  ( X )  • 

Squint  angle  (n0)  * 
Slant  range  (rQ)  “ 

The  half  aperture  length  is  given  by 


3max 


XrQ 

4psinnp 


1.358  m, 
0.0321  m  , 

6*. 

100  km . 


With  these  parameters,  the  quadratic,  cubic  and  quartic  phase  terms  at  the 
end  of  the  full  aperture  are  given  by 


quadratic  phase 


4  it  2 

2smax' 


cubic  phase 


quartic  phase 


irXrQ 


■  3 1 . 1  tt  rad. 


4tt 

3 

—  a 

X 

3smax» 

ir  A 

2rF 

32p  3 

tannp 

12. 

2  it  rad 

4tt  4 
*  ~d43max • 
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n\3rp 

51204 


‘  ?. 

tan  np 


-o. 


■  0.686n  rad. 


In  practice,  the  overall  phase  error,  with  a  heavil/  weighted  aperture, 
should  not  be  more  than  0.75tt  rad.  It  can  be  seen  then,  that  even  in  this 
extreme  case,  the  contribution  due  to  the  quartic  phase  term  can  be 
ignored. 
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APPENDIX  B 


Azimnth  Resolution 


The  azimuth  resolution  is  mainly  determined  by  the  second  order 
phase  term  given  in  (3.2). 

Let 


9 


a  ;_3 


2 


The  Doppler  bandwidth  (BW)  is  given  by 


1  d9 

BW  * - 

2ir  ds 


3-NpASA 


l2NpASA, 


2sinnp 

Ai*j? 


M  A  C 
"P“WA ' 


where  all  variables  are  defined  in  Section  3 . 

The  resolution,  p',  in  the  flight  direction  is  given  by 
1.4  1 .4  Arp 

p*3 - - 

BW  2NpASAsm  rip 


where  the  factor  of  1  .4  is  due  to  the  broadening  effect  of  the  Hamming 
window. 


The  azimuth  resolution  in  the  direction  perpendicular  to  the 
antenna  pointing  vector  is 

1 .4  Arp 

p  =  p'sinnp  •  - . 

2NpASAsinnF 
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